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ABSTRACT 


The form and size of a channel in cohesionless material, 
stable against erosion for a definite discharge, Q, are studied. 
The angle of internal friction@ and the limiting tractive 


force t 4, are taken as known. Distribution of shearing stresses 


t is assumed to be such that they are proportional to the distance 
between bottom and water surface, measured at right angles to 
the bottom. In addition to the action of gravity and ry 
stress t the grains are acted upon by a hydrodynamic lift force, 
proving to be proportional to t. The differential equation of the 
bottom form is established and integrated numerically; the form 
depends on@. 

Based on the logarithmic law of velocity distribution and 
the assumed distribution of shearing stresses, the velocities in 
all parts of the cross section can be found, and the total dis- 
charge is found by numerical integration. ‘ 

A profile consisting of the curved bank-part of the above 
mentioned cross section and a“middle-part" of indefinite width and 
of constant depth y would be stable for the same tractive 
force. On the assuffiption however that nature will produce that 
cross section which has a minimum of area, only one definite 
solution, viz. the equilibrium profile, is found. The dimensions 
depend not on@ alone but also on the relative roughness of the 


bottom . Provided that the hydraulic roughness k is assumed 


max 
to be in conformity with that of natural watercourses, it is found 
that the area of the equilibrium profile varies slowly with a and 
must be proportional to 
0.9 
max 

The above assumptions are checked by calculation of a 
complete set of isovels. 

Further three model tests, carried out in Vienna in 1916, 
are studied and compared with profiles calculated according to 
this theory. The values of 7p are found to be varying from 14° to 


20": On the same basis a study is finally made of the relation 
v 
between mean and maximum velocities, =, resulting in a 
max 


simple diagram giving — as a function of 0, the “degree of 
Ymax 


fullness” of the profile, and also as a function of ——» the 


reciprocal relative roughness. 
Methods for estimating k are given. 
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STABLE CHANNELS 


A. E. Bretting 
Professor of Hydraulics, 
Technical University of Denmark, 
Copenhagen, Denmark. 


1, INTRODUCTION 


The term “stable channels’ has generally been used as 
relating to the cross section of a channel which will be stable 
for a definite discharge both against erosion and against sediment- 
ation. In this paper only the question of stability against erosion 


has been treated, and only as far as cohesionless bottom material 
is concerned, 


The subject of this investigation is the cross section ultimately 
created when nature itself, at a definite discharge Q (m* /sec), 
excavates a channel bed in sand material with a definite limiting 
tractive force t (kg/m*) just able to move the sand grains on 
a horizontal botf@f and with an angle of internal friction of so 
degrees. 


The specific weight of the water is y (kg/m*), and the 
equivalent sand-roughness of the bottom of the channel is k (m); 
this roughness is taken to be uniform over the whole width of the 
channel, 


It is further assumed that the velocity v_ at every point P of 
the cross section can be found according to thé usual logarithmic 
law of velocity distribution 


v 
Zz 
= 8.48 + 2.5 In 
where k is the roughness of the bottom (m), 


z is the distance from the bottom (m) measured in the 
direction perpendicular to the bottom, 


1 
| 
| 


2 SEC. 1. 


st (m/sec) , 


= the shearing stress (kg/m*) at the bottom in a line 
through the point P perpendicular to the bottom, 


9 = specific weight of the water (kg/m* ) 


g, acceleration of gravity = 9.81 m/sec’. 


To find the velocities, a definite law for the distribution of the 
shearing stresses along the bottom contour must be assumed. 


In the middle part of the cross section, where the bottom line 
is horizontal, the shearing stress will have its maximum value 
t equal to the limiting tractive force of the sand grains in 
qit#tion. 


This limiting tractive force will strictly speaking depend on the 
longitudinal slope of the channel. 


The relation between the limiting tractive forces t_ for a plane 
bottom sloping in the direction of the flow at an angle v and Th for 
a horizontal bottom, can be shown to be: 

t 
= cos v(l - tg v cot), 


where @ is the angle of internal friction of the sand. 


Since the longitudinal slope = sin v will generally be small we 
take the limiting tractive force to be independent of the slope. 


The following law for the variation of + is assumed: 


’ 


max y max’ 


where y is the depth of water at the point in question, and & is the 
inclination of the bottom with a horizontal plane. (Fig.1). The 
shearing stress will consequently be proportional to the length of a 
normal to the bottom reckoned between the bottom and the water 
surface. The shearing stress on the slopes will be somewhat greater 
than that found by using the hypothesis that + varies proportionaliy 
with the depth y , an assumption previously used. Transfer of 
shearing forces from the middle of the section against the banks is 
thus to a certain degree taken into account. 


A justification of this assumption is later found by calculation 
of the corresponding complete set of isovels (Fig.7 ), which are in fair 
accordance with experience from actual measurements. 


i Vv — 
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2. ESTABLISHMENT OF THE DIFFERENTIAL EQUATION 
FOR THE SHAPE OF THE PROFILE. 


2.1. EQUATION ON THE ASSUMPTION L = c- 1. 


The profile sought should just be in equilibrium at every point 
of the bottom, given the assumed distribution of shearing stresses. 


On the bottom an element of one unit area will be stressed 
by the force 1 in the direction of the flow. The longitudinal slope of 
the channel is considered insignificant. The surface element has the 
inclination 2° with a horizontal plane. 


The submerged weight of sand grains per unit area stressed 
by +t is designated as W. (Fig.2). This force is acting vertically 
downwards and is resolved into the forces W-cos @ in the direction of 
the normal to the element, and W-sin © acting in the plane of the 
element in the direction of its transversal slope. 


The grains of the element are further acted upon by the hydro- 
dynamic lifting force L, which is upwards directed in the normal to 
the plane of the element. The resulting force in this direction will 
consequently be (W-cos & - L), acting downwards. 


The total stress on the element in question will be the resultant 
of the three above-mentioned forces : 
1) + in the plane of the element and in the direction of the flow, 
2) (W-cos 4- L) in the direction of the normal to the element, 
3) W-sin @ in the plane of the element perpendicular to the direction 
of flow. 


The resultant of forces 1) and 3) is\x? + W* . sin’a acting 
in the plane of the element, whereas force 2) is perpendicular to the 
said resultant. 


If the angle of friction of the sand is taken to be @ degrees, 
it is a condition of equilibrium for the sand grains of this element 
that + - sin’ a 
(2) ‘cosa - L 


The magnitude of the hydrodynamic uplift L will be studied 
below. 


We put provisionally 
(3) Lee-t, 
where c will prove to be a constant (Section 2.2). 
From equation (2) we get : 
(2a) W*+cos* a - 2-L»W-cosa + = + W*: sin® a) 


| 
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4 SEC. 2.1, 


From equation 1,(1) we get :; 


t 


copa 
Fory = have a = 0° , cosa = 1 Bina = 0 
Tmax ’ Which, inserted in (2a), gives 
W -2crt 


max W +c cot? ; taking the 
square root of both sides we obtain : 


(4) Ws (c + cot@P) 


This value of W is inserted in (2a) together with L from (3), 
which gives : 


(5) 2 cos a 


+ cot?® )-cos* a - cot” P| = 0 
Substituting 1, (1); 


cos a 


we get from (5) : 


+ + cot’p )-cos® a - } = 0 
The following substitutions are used : 
d d 
(6a) =n 3 (6b) = (6c) tan = = 3 
(6d) — ey 


cos* a 


(6a) - (6d) are inserted in (5a) : 


1+ 


3 
max max 
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- (1 + + cot*p -(1 + 


+ (1 + + SRP) 


This is the correct differential equation in 7 and § for the 
profile in question. It can be written in the form 


4 

c c [ 
cotte (1+ SEP) + - 


2.2 DETERMINATION OF THE FACTOR c IN FORMULA 2,1 (3) 


In the paper by H. A. Einstein (Ref. 1] information is found 
bearing upon this subject. The notations of Einstein’ s paper are 
used ( in this section only ). Formula (36), page 31 runs 

2 


Py = Cy; 5 where 


Py, = average lift pressure per unit of area, 
0,178 

By = density of the fluid = oP 

¥, = specific weight of fluid, 

g = acceleration of gravity , 


u__= the flow velocity at a distance 0. 35°Dye from the theoretical 
bed, 


Dys = sieve size of the grains of which 35 percent are finer. 


It is indicated that the pressure fluctuations due to turbulence 
in their duration follow the normal error law, the standard deviatian 
being 0. 364 of the average lift. 


A deviation from the mean of 2.75 times "standard deviation , 
viz. 


2.75 0.364 - = 1.0. 


PL, PL, 
will only have a statistical probability of 6 per thousand to be 
exceeded and corresponds to a practical maximum value of the 
hydrodynamic lift force per area : 

L= 2p ‘ which is very seldom exceeded. For greater 
values of eynélds’ numbers, which are exclusively considered, is 
found 


5.75- log ( 30.2 , where 


3 
) 
| 
3 
4, 

| 
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y = 0.35(0.77 +> k_), 
A: k the equialent sand roughness, 


u. where 


a 


shearing stress along the bottom 
usu, + 5.75 log (30.2 + 0.35 + 0,77) 
u = 5.24 


The mean value p; according to [1] formula (36) 


PL = 2.45; = 2.45° = 2.45 
f 


We consequently find the maximum value of the hydrodynamic 
lift force, defined as above, 


L=2p, = 


The constant c in equation (3) is found to be 
(9) c = 4.9 


3. SOLUTION OF THE DIFFERENTIAL EQUATION 2.1(8) 
In the numerical calculations of ¢ for varying values of », 
(a ) in the last te of the numerator can be 


omitted from the first approximation,whereupon the said last term 
is corrected. 


3.1. FORMULAE FOR NUMERICAL INTEGRATION 
The following substitutions are used : 


2.2.(9) c= 4.9 

(10) 

(11) b = cot?@+(1 + cote =a: 
(12) m= 1 - 


an | the differential equation (8) takes the following form : 


b+2n-2m- 


‘ 
a 
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For a given value of P the — a, b and m are fixed 


constants, and the numerical values of re can be found 
corresponding to chosen values of n. 


The integration of. the differential equation (13) is made 
numerically for values of @ = 15°, 20°, 25° and 30°, 


The corresponding values of the constants a, b and m are 
found from formulae (9), (10), (11) and (12), and the following 


Ter forms of (13) are found : 


dn ®  1-761643 -2 + 0.238357 - {1 + 
(13a) 7 
24,536516 + 2 - 2+ 0.238357 - 
dn 2 1.56071 - 2n+ 0.43929 - + 
(1 3b) (re 
11,78125 + 27-2. 0. . 
= 25°) | 
dn 1.43766 - 2 + 0.56234 fi + 
(13c) 
6.61168 + 2n- 2° 0.56234 


(3d) (GL) 


2 1.35348 - 0.64652 [i + 


4.06044 + 2n- 2° 0.64652 - 


As values of 4 are taken : 
0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 
0.75, 0.8, 0.85, 0.9, 0.95, 0.98, 1.0. 2 
For each value of @ and these values of » the values of Cc ) are 


calculated from (13a) b) c _and the corre rf ing values of 
d d 
d 


ables 3.1./la/b/c/d/ for p = 15°, 20°, 25° “and 30° resp. 


| 
if 
4 
30°| 
4 
in 


Table 3.1/1a/ 


a 

0 0.071797 0.267949 3.73205 1.035276 0.000000 
0.1 0.063240 0.251476 3.976529 1.031135 0.103114 
0.2 0.055030 0.234585 4.262854 1.027147 0.205430 
0.3 0.047149 0.217139 4.605354 1.023303 0.306991 
0.4 0.039582 0.198952 5.026360 1.019599 0.407840 
0.5 0.032312 0.179756 5.563085 1.016028 0.508014 
0.6 0.025328 0.159148 6.283474 1.012585 0.607551 
0.7 0.018616 0.136440 7,329257 1.009265 0.706485 
0.75 0.015358 0.123927 8.069238 1.007650 0.755737 
0.8 0.012164 0.110291 9.066962 1.006064 0.804851 
0.85 0.009032 0.095039 10.521951 1.004506 0.853830 
0.9 0.005962 0.077215 12.950913 1.002977 0.902679 
0.95 0.002952 0.054329 18.406283 1.001475 0.951401 
0.98 0.001174 0.034259 29.189274 1.000587 0.980575 
1.0 0.000000 0.000000 oo 1.000000 1.000000 


| 


Table 3,1/1b/ 


2 092)" 

0 0.13247 0.36396 2.74748 1.06418 0,000000 
0.1 0.114025 0.33768 2.96142 1.05547 0.105547 
0.2 0,097023 0,31149 3.21042 1.04739 0, 209478 
0.3 0.081328 0,28518 3.50655 1.03987 0, 311961 
0.4 0.066822 0.25850 3.86848 1.03287 0.413148 
0.5 0.053404 0,23109 4.32726 1.02635 0.513175 
0.6 0.040988 0,20245 4.93937 1.02029 0.612174 
0.7 0.029500 0.17176 5.82222 1.01464 0.710248 
0.75 0.024083 0.15519 6.44384 1.01197 0.758978 
0.8 0.018875 0.13739 7.27874 1.009393 0.807514 
0.85 0.013869 0.11777 8.49136 1.006910 0.855874 
0.9 0.009057 0.095168 10.50770 1.004518 0.904066 
0.95 0.004436 0.066606 15.01359 1,002216 0.952105 
0.98 0.001752 0.041862 23.888195 1.000876 0.980858 
1.0 0.000000 0,000000 oo 1,000000 1.000000 


! 
j 


10 


= 25°] 


Table 3,1/1c/ 


ag an 

0 0.21744 0.46631 2.14451 1.10338 0,000000 
0.1 0.182852 0.42761 2.33857 1.08759 0.108759 
0.2 0.152249 0.39019 2.56285 1.07343 0.214686 
0.3 0.125035 0.35360 2.82803 1.06068 0.318204 
0.4 0.100742 0.31740 3.15061 1.04916 
0.5 0.079002 0.28107 3.55780 1.03875 0.519375 
0.6 0.059514 0.24395 4.09912 1.02933 0.617598 
0.7 0.042047 0.20505 4.87677 1.02081 0.714567 
0.75 0.034008 0.18441 5.42262 1.01686 0.762645 
0.8 0.026404 0.16249 6.15411 1.01312 0.810496 
0.85 0.019216 0.13862 7.21387 1.00956 0.858126 
0.9 0.012429 0.11149 8.96978 1.00620 0.905580 
0.95 0.006027 0.07763 12.88099 1.00301 0. 952860 
0.98 0.002367 0.04865 20.55420 1.00118 0,981156 
1.0 0,000000 0.000000 o 1.00000 1.000000 


: 


= 30°, Table 3,1/1d/ 
an 2 2 

0 0.333333 0.577350 1.73205 1.154703 0.000000 
0.1 0.273203 0.522687 1.91319 1.128362 0.112836 
0.2 0.222426 0.471620 2.12035 1.105634 0.221127 
0.3 0.179031 0.423121 2.36339 1.08583 0.325749 
0.4 0.141608 0.376308 2.65740 1.06846 0.427384 
0.5 0.109141 0.330365 3.02696 1,053158 0.526579 
0.6 0.080864 0.284366 3.51659 1.039647 0.623788 
0.7 0.056200 0.237065 4.21825 1.027715 0.719401 
0.75 0.045082 0.212325 4.709756 1.022293 0.766720 
0.8 0.034710 0.186306 5.36750 1.017207 0.813766 
0.85 0.025046 0.158259 6.31874 1.012446 0, 860579 
0.9 0.016059 0.126724 7.89116 1.007998 0.907198 
0.95 0.007718 0.087851 11,38297 1.003851 0.953659 
0.98 0.003015 0.054905 18.21316 1.001506 0.981476 
1.0 0.000000 oo 1.000000 1.000000 


0.000000 
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3,2, CALCULATION OF THE INTEGRAL, ESPECIALLY 


1-4 
Calculation of the integral 8 = fi n) + dn ie made by the 


0 
trapezoidal rule except for the last interval 0.98 4 31 , where 
this rule is not usable because f( 1) tends to infinity for y = 1. 
4 ‘ 


In this region ( g ) in the formula 3.1.(13) can be completely 
neglected, and we get 


(14) d§ \/b+2y-2m-: 


a-2n+m-: 
We substitute : 
nz=1-4 , =1- 24 +4" , where J is small in relation 
and find 


+ 2(l-m) -24(1-2m) - 2m -4’ 
(a- 2+ m) + 24(1-m)+ 
According to 3.1.(10) and 3,1,(12): 


a-2+ m=0_, consequently 
+ -m -m 
15 


In (15) the latter square root is expanded into a series, and 
considering that 4 is small compared with 1 we get, neglecting 
terms with 4 * and higher powers : 


: 

a 

-1/2 1. yg fm 44. __1-2m 

n -m I-m b+2(I-m) 

m 1-2m)? 

bF2(I-m) m(bt 2(m- 
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By integration we find ; 


m m- P 
The constants b and m are ner for each value of P from the 
above-mentioned calculations (compare the formula 3.1.(13) with 3.1. 


(13a), 3.1.(13b), 3.1,(13c) and 3.1, (13d). 

For , b= 24,536516 , m = 0,238357 we find : 

(16a) dn = - {1 - 0,032773.48 + 0. 00028184 
1-4 


For Z4= 0.02 


(16aa) as 1.169876-{ 1 - 0,00065546 + 0000000112 } 


0,98 
= 1.169109, 


For [p = 20> 49 = 11,78125 , m = 0,43929 we find 


(16b) dy = 6.78400°Va - 0.068427A + 
For A= 0.02 


(16bb) $3. 0, 950400- 1 - 0.00136854 + 0. 00000202} 


.98 
= 0.958089, 


For fp = 25°] b = 7.48700 , m = 0.56234 we find 
(16c) dn= 5, 84926-YA - 0.1015234 + 0.014835. 
1-A 


y 

d d § 
" an 

4 . 

1 

a 

0 

2 
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For & = 0.02 


1 
| dn* 0.827208: 1 — 0,00203046 + 0. 00000593} 


0.98 
= 0.825533. 


= 4.06044 , m = 0.64652 we find 


- 0.13194 + 0.02964 *} 


For Z = 0.02 


(16dd) dn = 0,734490° { - 0,02638 + 0. 00001184} 


. 98 
= 0.732561. 


3.3, NUMERICAL INTEGRATION. COMPUTATION OF n,§ 
VALUES FOR @ = 15°, 20°, 25° and 30°, 


The shape of the bank-part of the bottom, i,e. between 
n= 0 andy = 1, wherey = Ymax’ ©an now be computed ( Tables 
3.3. /2a/b)c/d/ for P = 15°, 20°, 25° and 30° respectively). 


For each interval in n the mean value of the two 
neighbouring values of d § is found, and 


n 

1{,d d d § 
For the last interval in n from 0.98 to 1.0 


a5 a$ dy, which are found in 
13,14 og * 


formulae3.2. (l16aa), (16bb), (16cc) and (16dd). 
The values of s n are found by successive addition of 4 s. 


: For |P = 30°|, > 

(16d) 
1-4 

| 
Ag 


Table 3.3/2a/ 


2 
0. 3.73205 0 
3.85429 0.385429 
0.1 3.97653 0.385429 
4.11969 0,411969 
0,2 4.26285 0. 797398 
4.43410 0.443410 
0.3 4.60535 1, 240808 
4.81586 0.481586 
0.4 5.02636 1, 722394 
5.29472 0.529472 
0.5 5. 56309 2.251866 
5.92328 0.592328 
0.6 6, 28347 2.844194 
6.80637 0.680637 
0.7 7.32926 3.524831 
7.69925 0.384963 
0,75 8.06924 3.909794 
8.56810 0,428405 
0.8 9, 06696 4.338199 
9.79445 0.489723 
0.85 10.52195 4.827922 
11,73643 0.586821 
0.9 12, 95091 5.414743 
15.67860 0, 783930 
0.95 18.40628 6.198673 
23.79777 0.713933 
0.98 29.18927 6.912606 
1.169109 
1.0 8.081715 


15 


5 
‘ 
A 
a 
| 


Table 3, 3/2b/ 


0. 2. 74748 0 
2.85445 0.285445 
0.1 2.96142 0, 285445 
3.08592 0.308592 
0.2 3.21042 0, 594037 
3.35849 0.335849 
0.3 3.50655 0.929886 
3.68752 0.368752 
0.4 3. 86848 1, 298638 
4.09787 0.409787 
0.5 4.32726 1. 708425 
4.63332 0.463332 
0.6 4.93937 2.171757 
5.38080 0.538080 
0.7 5.82222 2.709837 
6.13303 0, 306652 
0.75 6.44384 3.016489 
6.86129 0.343065 
0.8 7. 27874 3.359554 
7.88505 0, 394253 
0.85 8.49136 3.753807 
9.49953 0.474977 
0.9 10,50770 4.228784 
12.76065 0.638033 
0.95 15.01359 4: 866817 
19.45090 0.583527 
0.98 23.88820 5.450344 
0.958089 
1.0 B= 6.408433. 


| 

. 

| 

| 


lp = 25°, 


Table 3, 3/2c/ 


A? a 2/m 
0 2.14451 0 
2.24154 0, 224154 
0,1 2.33857 0.224154 
2.45071 0, 245071 
0.2 2.56285 0.469225 
2.69544 0.269544 
0,3 2.82803 0.738769 
2.98932 0, 298932 
0.4 3.15061 1.037701 
3.35421 0.335421 
0.5 3.55780 1,373122 
3.82846 0, 382846 
0.6 4.09912 1, 755968 
4.48795 0.448795 
0,7 4.87677 2.204763 
5.14970 6.257485 
0.75 5.42262 2.462248 
5.78837 0, 289419 
0.8 6.15411 2.751667 
6.68399 0, 334200 
0.85 7,21387 3.085867 
8.09183 0.404592 
0.9 8.96978 3.490459 
10.92539 0.546270 
0.95 12, 88099 4.036729 
16.71760 0.501528 
0.98 20,55420 4.538257 
0. 825533 
363790 


i 
17 
i 
| 
a 
- 
. 
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IP = 30°. | Table3, 3/2d/ 
a (<3) 

an Im 3 

0 1, 73205 0 
1.82262 0.182262 

0.1 1.91319 0.182262 
2.01677 0.201677 

0.2 2.12035 0.383939 
2.24187 0.224187 

0.3 2.36339 0. 608126 
2.51040 0.251040 

0.4 2.65740 0. 859166 
2.84218 0.284218 

0.5 3.02696 1.143384 
3.27178 0.327178 

0.6 3.51659 1.470562 
3.86742 0.386742 

0.7 4.21825 1.857304 
4.46401 0, 223200 

0.75 4, 70976 2.080504 
5.03863 0.251932 

0.8 5.36750 2.332436 
5.84312 0.292156 

0.85 6.31874 2.624592 
7.10495 0.355248 

0.9 7,89116 2.979840 
9.63707 0.481854 

0.95 11,38297 3.461694 
14.79807 0.443942 

0.98 18,21316 3.905636 

0.732561 
1.0 = 4,639197 | 
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4. DETERMINATION OF PERIMETER, AREA 
AND DISCHARGE FOR "BANK-PART" AND 
“MIDDLE-PART. PRINCIPLE OF MINIMUM 
WORK. 


4.1. INTRODUCTION, 


The profile determined in the foregoing section should 
be stable against erosion for every value of » between zero and one 
as long as the limiting tractive force is not exceeded. The relative 
width of the curved sloping bottom at one side of the axis of the 
channel is called 6 = 1.0 


The absolute width of this part of the cross section is 
B.Yy , and the corresponding parts of the perimeter, area and 
discha¥e are called Up: F, and Qp, respectively. This part of the 


cross section is called the “bank-part’ . (Fig.3). The shearing 
stress in this part varies between + = 0 for($ = 0, » = 0) and 
T= T for ( § = 6B, n= 1), and it is evident that another half 
cross 88&etion consisting of the said “bank-part” and a “middle- 
part" with constant depth Ymax 2nd an arbitrary width Bo.y. 4, 


between the " bank-part" and the axis should be equally stable. 
The perimeter, area and discharge of such a “middle-part* are 
called Uso: Fao and respectively. 


So far an infinity of possible solutions seem to exist for 
fixed values of the total discharge of the half-section Q = Q, 7 Qe 


and for and @. 


If we imagine, however, that the profile is eroded 
gradually in a uniform mass of sand at a constant discharge Q, the 
intensity of the erosion will steadily diminish, and the profile must 
asymptotically approach the state of equilibrium commensurate 
with the sand grains in all parts of the bottom being in a state of 
incipient motion, 


It thus seems natural to assume the principle of minimum 
of work and assume that nature will produce the profile that 
requires a minimum of erosion, i.e. a minimum of the cross 
section F = F, + F According to this principle only one definite 


equilibrium profile exist for fixed values of Q, and @. 


We therefore proceed to find the cross sections Fs and 
F and the discharges and where_upon the indeterminate 
' quantity By is fixed in such a way that it causes F = Fy + Feo 
to be a minimum for a fixed value of Q = Q, + Qe or the total 


discharge in the half cross section, and for fixed values of 1 
and @ . max 


: 

4 

‘ea 

oe 

3 

4 

id 

4 

if 
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4.2, DERIVATION OF FORMULAE FOR MEAN VELOCITY, AREA 
AND DISCHARGE OF AN ELEMENT OF THE CROSS SECTION, 


As previously mentioned, it is assumed that the loga- 
rithmic law of velocity distribution is valid, and that the velocity 
at a point P, situated at the distance z from the bottom, and 
measured perpendicularly to the bottom element ds, will depend 


exclusively on the shearing stress + on this element and on the 
distance z, 


The length of the normal to the water surface is 


Fig. 4) 
From equation 1,(1) we get by means of (17) : 
t y 29 


The values of § are compiled in Tables 3,1/1a/b/c/d/ 
We consequently have ; 


(19) 


max max 


z 
= 8.48 + 2.5 
where k is the equivalent sand roughness of the bottom supposed to 
be constant for the whole width of the channel. 


If the element of the cross section shown in Fig.4 had 
a constant width equal to the base ds, the mean velocity in such 
an element could be found as the velocity at a point 0.3679 - Zo 


from the bottom. The curvature of the bottom is slight, so that it 
is considered permissible to disregard this and find, at the same 
distance, the mean velocity in the actually wedge-shaped element 
; but for the computation of the element area and its dis- 
Fye the wedge-shape must be taken into account. 


We consequently get : 


{0.48 + 2.5 In ( 0.3679, 


2.5V/ VE {in + In 


The area of the surface-element (Fig.4), with base ds, 


| 
| 
| 
v 
[(20)] 
=> 
mes 
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measured between the bottom and the water-surface, and when r 
is the radius of curvature of the bottom, will be 

(21) dF = + de-(1 + ) = de(1- 


For the radius of curvature r we have 


or from 2,1, (6a)b)c)d) 


(22) 


y -a§ 2 


(17) Ymax 


From (21), (22), (23) and (17) we finally get 
(Sa) 


(24) =: (a$yyi + + . 


2 
max 1 + ( day 
By the following calculations with bray differences we 
introduce the mean values of § and , * °&§ | for the sides 
of the element in question and put 4 ford . and 4 a) for 


a (24) then takes the form 


| 

1+ (x) 

x 

d? 

1 + ( ) 

372 

4 

| 
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where subscript m signifies mean values. 


Accordingly the mean velocity in the element considered 
is taken from (20) 


The discharge 4 Q in the element 4 F is; 


2.5 Y max 


4.3. FORMULAE FOR THE QUANTITIES S AND FOR DIMENSION- 
LESS VALUES OF PERIMETER, AREA AND DISCHARGE U’, 
¥: AND QQ’ RESPECTIVELY, VALID FOR THE “ BANK- PART. 

"MIDDLE-PART" AND FOR THE TOTAL HALF-SECTION. 


The following notations are introduced : 


(26) Sq * 
1 
2 
(27) 1 + Gy 
n=0 
nel 
2 y2 


d 
a 2(1 + Cr 
m 
7=0 


(29) 


(30) 


(31) 


(32) 


(33) 
(34) 


(35) 
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‘rg? 
1 
alse 
m 
(42) 
53,9 yop in 


So = - 
Ss = 83.1 - 53.2 


For the “bank- part” of the cross section we find the 


following dimensionless expressions for perimeter, area and 
discharge. 


23 

n=O 2, 1+ 

m 

U 
0 

Ymax 

y 

max 

11 

~ (Sp, 1 8p, 2) Int 

‘ 

2.5 9 nax 

+ 

| $3.1 53,2 

11 

: 


24 


(26a) 
(27a) 
(28a) 
(29a) 


38a 
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For the “middle-part” of the section we have 


21 


nd thereby 
So Bo 
= Bo (30a) Si 2 = 0 
8,70 (32ay S, = 0 

= Bo 
Ymax 

2 0 
max 

Q lly 

B, 0 max 


max 2 
2.5 * Ymax 
We put for the total half-section of the channel 


U= U, + Ung 


F = F, + Feo 


Q = Q, + Q, o, and find from (36), (37), (38), (36a) 
(37a), (38a) 


U 
» 0 0 
max 
Q + B in(7™8X)4 
= (Sp + Bg) 83. 
SS 


q 
| 
| 
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4.4. DETERMINATION OF 68, CORRESPONDING TO MINIMUM 
OF CROSS SECTION, 


For a fixed value of @ the quantities S;, S, and S5 


are constants, and if further Q, Vides 9 and k are constant, 
equation 4,3.(44) represents a relation between Bo and y as 
which must be satisfied : ” 
lly 
Q max 
(44a) we— - (S, + By) In( ——)- S, = 0 
3 
max 
together with 
(43a) F = (8, + Bp) Ymay 
F is a function of 6, and | ee and to obtain a minimum 
value of F simultaneously with wy = 0, Layrange’s method 


gives the following condition ; 


0 = °Ymax hence 

2(S; Bo) max 


- (S, + Bg) Inf + (Sp + By) = 0. 
2.5 max 2 


When deducting (44a) the result is: 
1 
By = 2(8, - Sp) In(—-™**)- (s, + , 


which gives the value of Bo corresponding to a minimum of F, 


2 3 max 

gy 

‘ 
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4.5. FORMULAE VALID FOR THE OPTIMAL PROFILE FOR 
DIMENSIONLESS VALUES OF PERIMETER, AREA, DIS- 
CHARGE, MEAN VELOCITY, HYDRAULIC RADIUS, WIDTH 
OF WATER SURFACE, SLOPE, SHEARING VELOCITY AND 
FULLNESS OF THE PROFILE,ALL EXPRESSED BY THE 
QUANTITIES S AND 


It follows from the —e section 4,4, that 
(46) S, + By = 2(8, - - 2 and 


(47) 9 + By = 2(S; - in(— - 28,. 


, (46) and (47) are introduced in 4.4. (42), 
4.4.(43) and 4,4. (44) in: 


= 5 = (Sy - - 2 Ss) 


max 


F’ = 2(8, - So) in( (S, - 8, 2 Sg) 


Ymax 


(5) -— = 2(8; 
2 


ly, 
max 
Equations (48), (49) and (50) are valid for the profile 
with a minimum of cross sectional area. 


For the mean velocity 


Vin = fr We get, dividing (50) by (49) : 


m 1 
Ym (51) 
9 
2(S, - S (— 2s, AX), S, 


2(S, - So) (S, - Sy - 2 S3) 


We want to introduce the slope of the water surface I 
(uniform motion assumed) and the hydraulic radius R of the cross 
section and find that the total shearing force K for the half cross 
section F must be : 


i 3 
| 
| 
| 
v 
: 
id 


7 
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(52) yFi= tds; we further have 


4.2(23a) As -a§ (V1 + (2 * Ymax 


1 
dn}. 
tml 4s St | 
and from 4,3(27) 
(53) ds = 1,1 + Ymax* "max - 


As s* 98 (52) gives with By from 4,4. (45) 


FI 
[54)] 51,1* fo 48, - 


* 


For the hydraulic radius R we find by 4,3. (42) 
R’ = R = F = 

y Ury 2 Sp + 
max max (So + Fo) 0 0 


and by (49) and 4.4. (45) 


Sa) (8, - 8, -2 


For the shearing velocity v, we have 


So + By 


t—s—) 0 
Ti 
2(8, - (Sy - S, - 2 Ss) 


max 
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For the width B of the half profile at the water surface 
we find 

(57) B’ = = (B+ Bo) = 2(S, - S») in(—-——)+ (B-S,-2 83). 


We later apply the term 2€, "degree of fullness. of the 
cross section, which is defined by 
11 


Ymaxy, 
B B+ ** Ymax 
2(S,-S) In(—, +(-S, -28,) 


found by using (49) and (57). 
4.6. EVALUATION OF NUMERICAL VALUES OF THE 
QUANTITIES S. 


It remains to compute the_numerical values of 
Sp S,, Sg and S, for@ = 15°, 20°, 25° and 30°. 


The calculations are given in Tables 4.6. /3a/b/c/d/. The 


fundamental values of V1 + 0 and ( 1+ Gt) 


are to be found in Tables 3.2./la/b/c/d/, whereas & § is found 
in Tables 3.4./2a/b/c/d/. The calculations are made from the 
formulae 4.3.(26) through 4.3. (35). 


The results are compiled in Table 4.6. /3e/. 


| 
4 
a 


- 15°, Table 4.6/3a/ 
| (22)? Sm 
0 1.035276 0 
0.051557 0.385429 1.033206 
0.1 1.031135 0.103114 
0.154272 0.411969 1.029141 
0.2 1.027147 0.205430 
0.256211 0.443410 1.025225 
0.3 1.023303 0.306991 
0.357416 0.481586 1.021451 
0.4 1.019599 0.407840 
0.457927 0.529472 1.017814 
0.5 1.016028 0.508014 
0.557783 0.592328 1.014307 
0.6 1.012585 0.607551 
0.657018 0.680637 1.010925 
0.7 1.009265 0.706485 
0.731111 0.384963 1.008458 
0.75 1.007650 0.755737 
0.780294 0.428405 1.006857 
0.8 1.006064 0.804851 
0.829341 0.489723 1.005285 
0.85 1.004506 0.853830 
0.878255 0.586821 1.003742 
0.9 1.002977 0.902679 
0.927040 0.783930 1.002226 
0.95 1.001475 0.951401 
0.965988 0.713933 1.001031 
0.98 1.000587 0.980575 
0.990288 1.169109 1.000294 
1.00 1.000000 1.000000 


To be continued 


fi 


Table 4. 6/3a/ (continued) 
%, 

0 

0.398228 0.020531 0.227062 0.004662 2.96507 
429074 0.025600 1.86004 
0.454505 0.116472 0.806173 0.058055 1.36176 
9.401917 0.175819 0.597843 0.105112 1.02886 
9.598004 0.246779 0.676703 0.166996 0. 78104 
0.335117 0.746849 0.250282 0, 58379 
9.688078 0.452076 0.810566 0.368437 0.42004 
0.283831 0.855050 0.242690 0.31319 
gsi 0.431343 0.336574 0.883343 0.297310 0.24808 
0.910682 0.371826 0.18713 
» ©, 817908 0.937153 0.484796 0.12982 
165675. 0.728382 0.962829 0.701278 0.07576 
yee 0.714669 0.690362 0.982847 0.678520 0.03461 


To be continued 
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31 

p= 15 Table 4.6/3a/ (continued) 

22 dn an)? 
0 0.267949 1.071797 

0.013823 0.016473 2.135037 
0.1 0.251476 1,063240 

0.048016 0.016891 2.118270 
0.2 0.234585 1.055030 

0.080283 0.017446 2.102179 
0.3 0.217139 1.047149 

0.108146 0.018187 2.086731 
0.4 0.198952 1.039582 

0.130431 0.019196 2.071894 
0.5 0.179756 1.032312 

0.146112 0.020608 2.057640 
0.6 0.159148 1.025328 

0.153918 0.022708 2.043944 
0.7 0.136440 1.018616 

0.076008 0.012513 2.033974 
0.75 0.123927 1.015358 

0.073757 0.013636 2.027522 
0.8 0.110291 1,012164 

0.069580 0.015252 2.021196 
0.85 0.095039 1.009032 

0.062936 0.017824 2.014994 
0.9 0.077215 1.005962 

0.053129 0.022886 2.008914 
0.95 0.054329 1.002952 

0.023484 0.020070 2.004126 
0.98 0.034259 1.001174 

0.011248 0.034259 2.001174 
1,00 0.000000 1.000000 


To be continued 


1, 050871 
= § 
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= 15°, Table 4. 6/3a/ (continued) 

0 

0.007716 0.002658 0.000021 0.000005 0.000015 
0.001113 0.000665 0.000684 
0, 200007 0.001943 0.001315 0.001027 
211898 0.003116 0.002327 0.001358 
wae 0.431673 0.004796 0.003887 0.001633 
0.608859 0.004095 0.003617 0.000897 
0. 0.004726 0.000884 
788 0.771332 0.006823 0.006394 0.000830 
9011892 0.889408 0.009700 0.009428 0.000714 
9-95 9.010014 0.983183 0.000344 0.009184 0.000818 
9.017119 0.980670 0.016788 0.016706 0.000168 
$0. 087042 40. 061414 

° 55,3 


p 
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5 


Sm 


0.75 
0.85 


0.95 
0.98 


1, 06418 
1.05547 
1.04739 
1.03987 
1.03287 
1.02635 
1.02029 
1.01464 
1.01197 
1.009393 
1.006910 
1.004518 
1.002216 
1.000876 


1.000000 


0 


0. 


0 


105547 


. 209478 
. 311961 
413148 
. 513175 
. 612174 
. 710248 
. 758978 
. 807514 
. 855874 
. 904066 
. 952105 
. 980858 
. 000000 


0. 


0 


0. 


052774 


. 157513 


260720 


. 362555 
. 463162 
. 562675 
. 661211 
. 734613 
. 783246 
. 831694 
. 879970 
. 928086 
. 966482 
. 990429 


.583527 1.001546 


- 285445 1.05983 
308592 1.05143 
335849 1.04363 
. 368752 1.03637 
.409787 1.02961 
-463332 1.02332 
.538080 1.01747 
. 306652 1.01331 
. 343065 1.01068 
. 394253 1.008152 
.474977 1.005714 
. 638033 1.003367 


.958089 1.000438 


To be continued 


a9 
0 | 
0,1 
0.2 
0.3 0 
0 0 
0.4 0 
0 0 | 
0.5 0 
0 0 
0.6 0 
0 0 
0.7 0 ; 
0 0 
0 
0 0 
0 0 ; 
0 7 
0 0 
0.9 0 
0 0 : 
0 
0 
0 
0 
j 
: 
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P = 20. Table 4.6/3b/ (continued) 
as’= 

2 45" Sin VS. As in§,,| 
0 

0.302523 0.015965 0.229726 0.003668 2.94173 
an 0.324463 0.051107 0.396879 0.020283 1.84824 
ie 0.350502 0.091383 0.510607 0.046661 1.34431 

0.382164 0.138555 0.602125 0.083427 1.01458 
0.421921 0.195418 0.680560 0.132994 0.76968 
sg 0.474137 0.266785 0.750117 0.200120 0.57505 
as 0.547480 0.362000 0.813149 0.294360 0.41368 
oi 0.310734 0.228269 0.857096 0.195648 0.30841 
0.346729 0.271574 0.885012 0.240346 0, 24430 
a 0.397467 0.330571 0.911973 0.301472 0.18429 
ie 1 0.477691 0.420354 0.938067 0.394320 0.12786 
0.640181 0.594143 0.963372 0.572381 0.07463 
nae 0.584429 0.564840 0.983098 0.555293 0.03408 
oie 0.958509 0.949335 0.995203 0.944781 0.00961 
1.00 

+ 6.518930 +4.480299 +3.985754 


To be continued 


x 

a 


20°, 


Table 4. 6/3b/ 


(continued) 
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To be continued 


0 0.36396 1.13247 

0.010790 0. 02628 2. 246495 
0.1 0.33768 1.114025 

0.037488 0.02619 2.211048 
0.2 0.31149 1.097023 

0.062727 0.02631 2.178351 
0.3 0.28518 1, 081328 

0.084643 0. 02668 2.148150 
0.4 0. 25850 1, 066822 

0.102363 0.02741 2.120226 
0.5 0.23109 1.053404 

0.115079 0.02864 2.094392 
0.6 0.20245 1.040988 

0.121771 0, 03069 2.070488 
0.7 0.17176 1, 029500 

0.060340 0.01657 2.053583 
0.75 0.15510 1.024083 

0.058717 0.01780 2.042958 
0.8 0.13739 1.018875 

0.055558 0.01962 2.032744 
0.85 0.11777 1.013869 

0.050418 0.022602 2. 022926 
0.9 0.095168 1.009057 

0.042717 0.028562 2.013493 
0.95 0.066606 1.004436 

0.018924 0.024744 2.006188 
0.98 0. 041862 1.001752 

0. 009079 0.041862 2.001752 
1.00 0.000000 1, 000000 

-0. 830614 

= 831 


= 

= 

4 

4 
wz 

ia 

a 

| 

j 

= 

| 

ite 
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Table 4.6/3b/ 


(continued) 


2 2 
Sin 


0.011698 
0.011845 
0.012078 
0.012420 
0.012928 
0.013675 
0.014823 
0.008069 
0.008713 
0.009652 
0.011173 
0.014185 
0.012334 
0.020913 


0.002785 
0.024810 
0.067975 
0.131446 
0.214519 
0.316603 
0.437200 
0, 539656 
0.613474 
0.691715 
0.774347 
0.861344 


0.000033 
0,000294 
0,000821 
0.001633 
002773 
. 004330 
006481 
. 004354 
005345 
. 0066 76 
. 008651 


012218 


0.934087 0.011521 


0.000008 
0.000117 
0.000419 
0.000983 
0.001887 
0.003248 
0.005270 
0.003732 
0.004730 
0.006088 
0.008115 
0.011770 
0.011326 


0.000024 
0.000216 
0.000563 
0.000997 
0.001452 
0.001868 
0.002180 
0.001151 
0.001156 
0.001122 
0.001038 
0.000878 
0, 000386 


0.980950 0.020515 0.020417 0.000196 


= Si 2 


+0.085645 +0.078110 -0.013227 


2 
0 
0.1 
: 0.2 
0.3 
0.4 
3 0 . 5 
0.6 
: 0.7 
0.75 
0.8 
0.85 
0.9 
0.95 
0.98 
4 
1,00 
2,2 3,2 | 
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2 


$, 


0 1,10338 0 

0.054380 0.224154 1.09549 
0.1 1.08759 0.108759 

0.161723 0.245071 1.08051 
0.2 1.07343 0.214686 

0.266445 0.269544 1.06706 
0.3 1.06068 0.318204 

0.368934 0.298932 1.05492 
0.4 1,04916 0.419664 

0.469520 0.335421 1.04396 
0.5 1.03875 0.519375 

0.568487 0.382846 1.03404 
0.6 1.02933 0.617598 

0.666083 0.448795 1.02507 
0.7 1.02081 0, 714567 

0.738606 0.257485 1.01884 
0.75 1.01686 0, 762645 

0.786571 0.289419 1.01499 
0.8 1,01312 0.810496 

0.834311 0.334200 1.01134 
0.85 1.00956 0.858126 

0.881853 0.404592 1.00788 
0.9 1.00620 0.905580 

0.929220 0.546270 1.00461 
0.95 1.00301 0.952860 

0.967008 0.501528 1.00210 
0.98 1.00118 0.981156 

0.990578 0.825533 1.00059 
1.00 1,00000 1.000000 


To be continued 


| 
| 
4 
4 
| 
i 
4 


Table 4.6/3c/ (continued) 
4s'= 

0 

0.245558 0.013353 0.233195 0.003113 2.91180 
a 0.264802 0.042825 0.402148 0.017222 1.82187 
~e 0.287620 0.076635 0.516183 0.039558 1.32259 
~ 0.315349 0.116343 0.607399 0.070667 0.99714 
“es 0.350166 0.164410 0.685215 0.112656 0.75604 
~e 0.395878 0.225051 0.753981 0.169684 0.56477 
“7 0.460046 0.306429 0.816139 0.250089 0.40634 
"e 0.262336 0.193763 0.859422 0.166524 0.30299 
my 0.293757 0.231061 0.886888 0.204925 0.24007 
“ 0.337990 0.281989 0.913406 0.257570 0.18115 
wi 0.407780 0.359602 0.939070 0.337691 0.12573 
Nii 0.548788 0.509945 0.963961 0.491567 0.07341 
ay 0.502581 0.486000 0.983366 0.477916 0.03355 
io 0.826020 0.818237 0.995278 0.814373 0.00947 

+5. 408671 +3.825643 +3. 413555 
= 8 8,1 


To be continued 
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Table 4:6/3c/ 
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(continued) 


2 

0 0.46631 1,217447 

0.009064 0.03870 2.400299 
0.1 0.42761 1.182852 

0.031376 0.03742 2.335101 
0,2 0.39019 1.152249 

0.052319 0.03659 2.277284 
0.3 0.35360 1,125035 

0.070465 0.03620 2.225777 

0.4 0.31740 1.100742 

0.085172 0.03633 2.179744 
0.5 0.28107 1.079002 

0.095832 0.03712 2.138516 
0.6 0.24395 1.059514 

0.101621 0.03890 2.101561 
0.7 0, 20505 1.042047 

0.050455 0.02064 2.076055 
0.75 0.18441 1.034008 

0.049196 0.02192 2.060412 
0.8 0.16249 1,026404 

0.046659 0.02387 2.045620 
0.85 0.13862 1,019216 

0.042458 0.02713 2.031645 
0.9 0.11149 1.012429 

0.036086 0.03386 2.018456 
0.95 0.07763 1.006027 

0.016034 0.02898 2.008394 
0.98 0.04865 1.002367 

0.007712 0.04865 2.002367 
1.00 0.00000 1,000000 

-0. 694449 

= 831 


To be continued 


‘ 
| | 
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Table 4.6/3c/ 


(continued) 


7 
2 


. 


0.016123 
0.016025 
0.016067 
0.016264 
0.016667 
0.017358 
0.018510 


0.009942 


0.010639 
0.011669 
0.013354 
0.016775 
0.014429 
0.024296 


0.002957 
0.026154 
0.070993 
0.136112 
0.220449 
0.323177 
0.443667 
0.545539 
0.618694 
0.696075 
0.777665 
0.863450 
0.935104 
0.981245 


0.000048 
0.000419 
0.001141 
0.002214 
0.003674 
0.005610 
0.008212 
0.005424 
0.006582 
0.008122 
0.010385 
0.014484 
0.013493 
0.023840 


0.000014 
0.000169 
0.000589 
0.001345 
0.002517 
0.004230 
0.006702 
0.004662 
0.005837 
0.007419 
0.009752 
0.013962 
0.013269 
0. 


023727 


0.000041 
0.000308 
0.000779 
0.001341 
0.001903 
0.002389 
0, 002723 
0.001413 
0.001401 
0.001344 
0.001226 
0.001025 
0.000445 
0.000225 


= § 


1,2 


+0.103648 +0.094194 


= Soe 


-0,016563 


= 83 9 


| 
| 

0 

‘ 0.1 

d 0.2 

0.4 

7 0.5 

4 0.7 

4 0.75 

0.8 

0.85 

0.9 

0.95 

4 0.98 

+4 
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1,15470 
1,12836 
1,10563 
1, 08583 
1, 06846 
1,053156 
1,039647 
1,027715 
1,022293 
1,017207 
1.012446 
1,007998 
1,003851 
1,001506 
1,000000 


0 

0.112836 
0.221127 
0.325749 
0.427384 
0.526579 
0.623788 
0.719401 
0.766720 
0, 813766 
0, 860579 
0.907198 
0.953659 
0.981476 
1.000000 


. 056418 
. 166982 
. 273438 
. 376567 
476982 


0 

0 

0 

0 

0 
0.575184 
0.671595 
0.743061 
0.790243 
0.837173 
0.883889 
0.930429 
0.967568 


0.990738 


0 
0 
0 
0 
0 
0 
0 


. 182262 
. 201677 
. 224187 
. 251040 
. 284218 
. 327178 
. 386742 
. 223200 
. 251932 
. 292156 
. 355248 
. 481854 
443942 
. 732561 


1,14153 

1,11700 

1,09573 

1.07715 

1, 060809 
1, 046403 
1.033681 
1,025004 
1.019750 
1.014827 
1,010222 
1.005925 
1, 002679 
1.000753 


To be continued 


| 
| 
0 
0.1 
0,2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.75 
0 
0, 8 . 
0 
0.85 
0 
0.9 
0 
0.95 
0 
0.98 
1,00 


ip = 30°, Table 4.6/3d/ (continued) 
as'= 

? VE, as’ 6% 
0 

0.208058 0.011738 0.237525 0.002788 2.87500 
0.23541 
0.00930 
1.00 

+4. 797025 +3, 374008 +3.018428 
= So 


To be continued 


| 
} 
} 


Table 4, 6/3d/ 


(continued) 


To be continued 


% ad dan 2 a 

0 0.577350 1, 333333 

0. 008016 0.054663 2.606535 
0.1 0, 522687 1, 273202 

0.027514 0.051067 2.495628 
0,2 0.471620 1, 222426 

0.045545 0.048499 2.401457 
0.3 0.423121 1,179031 

0.061028 0.046813 2. 320639 
0.4 0, 376308 1, 141608 

0.073525 0.045943 2.250749 
0.5 0.330365 1,109141 

0.082599 0.045999 2.190005 
0.6 0. 284366 1, 080864 

0.087592 0.047301 2.137064 
0,7 0.237065 1,056200 

0.043518 0.024740 2.101282 
0.75 0.212325 1.045082 

0, 042486 0.026019 2.079792 
0.8 0.186306 1,034710 

0.040364 0.028047 2.059756 
0.85 0.158259 1.025046 

0.036810 0.031535 2.041105 
0.9 0.126724 1.016059 

0.031369 0.038873 2.023777 
0.95 0.087851 1,007718 

0.013968 0.032946 2.010733 
0.98 0.054905 1,003015 

0,006723 0.054905 2.003015 
1,00 0.000000 1,000000 

-0.601058 


| 
is 
3 
4 
ad 
| 
< 
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ip = 30°, | Table 4.6/3d/ (continued) 

0 
0.020972 0.003183 0.000067 0.000016 0,000046 
a 0.020463 0.027883 0.000571 0.000233 0.000417 
0.020196 0.074768 0.001510 0.000790 0.001024 
ne) 0.020172 0.141803 0.002860 0.001755 0.001714 
waa 0.020412 0.227512 0.004644 0.003207 0.002374 
0.021004 0.330837 0.006949 0.005270 0.002915 
eT 0.022135 0.451040 0.009984 0.008182 0.003257 
aa 0.011774 0.552140 0.006501 0.005604 0.001664 
tT 0.012510 0.624484 0.007812 0.006945 0.001635 
ai 0.013617 0.700859 0.009544 0.008732 0.001552 
ee 0.015450 0.781260 0.012070 0.011348 0.001401 
wi 0.019208 0.865698 0.016628 0.016039 0.001157 
eT 0.016385 0.936188 0.015339 0.015088 0.000497 
we 0.027411 0.981562 0.026906 0.026781 0.000249 


+0.121385 +0.109990 - 0.019902 


= § 


2,2 


= 83 9 


| 
A 


Table 4.6/3e/ 


+8,081715 


45 


S14" +5. 535016 


Bs “8, -0.067042 
= +8.167180 S, = +5.467974 
+4.907011 8, -1.050871 
= +4,845597 8, = -1,040952 
+4.4802¥9 
B = +6. 408433 9" -0.085645 
S, = +6. 518930 S, = +4.394654 
+3.985754 S, -0.830614 
S, = +3,907644 S, = -0.817387 
B = +5, 363790 9° -0.103648 
Sy = +5. 498671 S, = +3. 721995 
+3.413555 S, -0.694449 
“Sy 9° -0.094194 +0.016563__ 
S. = +3,319361 S, = -0.677886 
30°, | S, +3.374008 
B = +4,638197 -0.121385 
S, = +4,797025 S, = +3. 252623 
+3.018428 S, ,= -0.601058 
“Sy = -0.109990 “Sy = +0.019902 
S, = +2.908438 S, = -0.581156 


‘ 
} 
| 
i 
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5. FORMULAE FOR Q’, F’, B’, R’ ANDI’ AS 


y 
FUNCTIONS OF p = -** For@ = 15°, 20°, 25° 
AND 30°. 


The fundamental formulae derived in section 4.5 for F’ 


and Q’ make it possible by simple division to find the formula 
forv’ . 


- Equally the formula for R’ can be found by division of 
the formulae for F’ and U’. 


We therefore give the formulae for Q’ and F’, R’ and 
I’ with the numerical values found for the quantities S and #6 for 
the different values of P . 


We further need the formulae for B’= (6 + Bp). 
From these the formula for Vy, can be found as 


v’ F 

« Vr. R 
and the formula for @ as 

ae = e 


The formulae are as follows : 


max 


Pp = 
(59a) Q’ = = 1,244754 1n?(11 p) 
2 max P y 2 
max 


+ 2.081904 In(1l p) - 1.040952 


(60a) F’ = 


F = = 1.244754 In(11 p) + 2.704281 
Ymax 
(61a) = B + By = 1.244754 In(11 p) + 5.318022 


R 1, 244754 ln(11 + 2.704281 
(S3a) n(Il p) + 5. 


(63a) yp = 1244754 1n(11 p) + 2.771323 
max 


Ymax * ( 9 


| 
| 
| 
: 
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(59b) QQ’ = = 0.974020 in*(11 p) 
+ 1,634774 In(11 p) - 0.817387 
(60b) = —E_ 0, 974020 in(11 p) + 2.121784 


(61b) B’ = + = 0.974020 In(11 p) + 4.135563 


R 0.974020 In(1l p) + 2.121784 


(63b) —- 9 974020 in(11 p) + 2.207429 . 


Ymax 
(59c) Q’ = — = 0, 805268 In’ (11 p) 
2.5 mea. y a 
max 


+ 1,355772 ln(1l p) - 0.677886 


(60c) . = 0,805268 In(11 p) + 1.758406 


max 
(61c) B’ = B+ By= 0,805968 p) + 3.400201 
R 0. 805268 1n(11 p) + 1. 758406 
62 
(63c) = 0.805268 In(11 p) + 1.862054. 


Ymax 


| 
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Q 


(59d) Q’ = 0.688370 In®(11 p) 
max 


+ 1.162312 In(11 p) - 0, 581156 


(60d) F’ « 


0.688370 p) + 1, 506497 


(61d) B’ = B+ By = 0.688370 ln(1l p) + 2.891576 


, R 0.688370 In(11 + 1,506497 
62d = 


(63d) = 0.688370 In(11 p) + 1.627882, 
max 
Numerical values of Q’, F’, (f+ Bo), R’ andl’ for 
fixed values of y 
max 
will be given later in the tables of section 7.1, viz. Tables 
7.1. /4a/b/c/d/ through 7.1. /8a/b/c/d/. 


6. DERIVATION OF FORMULAE OF THE 
pBe. 
TYPE z’ = Ag: p A,» WHERE 


max 
- AND Ag AND Be ARE 


FUNCTIONS OF @® ALONE, WHEREAS 
Ay IS FUNCTION OF p ALONE, 


6.1, INTRODUCTION, 


In formulae 5(59a)b)c)d) through 5(63a)b)c)d) the dimen- 
sionless quantities Q’, F’, B’, R’ and I’ are expressed as 


y 
functions of p = 


When z’ in general stands for any of these quantities 
we try to adjust formulae of the type 


to the numerical values of z’, which are computed for fixed 
values of p by the above-mentioned formulae for each value of - ‘ 


=" the numerical coefficients varying with Pp. 
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Putting in the first approximation A, = 1, we find for 


each value of P a simple power formula giving certain numerical 
values of A@ and Be . 


We then compute 
(65) log A, = logz’ - (log Ag + Boe logp), 
giving the errors of the power formula verying with p. 


It will then appear that these values of log A. are 
practically independent of @ . P 


Consequently we are now able to expressAg and Be 
as functions of @ alone and A, as function of p alone, 
6.2. (A, = 1). SIMPLE POWER FORMULA. 

For the first approximation with A= 1 we tabulate for 


each value of @ the corresponding values of z’ for the fixed values: 
p= 10, 20, 30, 50, 100, 200, 300, 500 and 1000 and take logarithms 


in (64) 
(66) log z’ = log Ag + Bo * log p. 


This form will on double logarithmic paper represent a 
straight line (Fig. 5). 


Its direction is taken parallel to the chord through the 
points 250 and Zo00' and it is placed so that the errors ey for 


p = 100 and p = 500 will have the same absolute value but 
opposite signs, which in the logarithmic representation corres- 
ponds to the same relative deviations for p = 100 and p = 500. 


It will prove that the devation for p = 24 will be nearly 
the same as for p = 500, and that this maximum deviation will 
not be exceeded between these limits. 


We find: = 


log 2; - log 2; 
200 50 _ 
Bg = fog = 1, 660964 (log 2599 - log 
and 1 
(68) log Ag = x (log 2500 + log 00) - 3.902411 


The numerical calculations are given later in the tables 
of section 7.1. 


6.3. DERIVATION OF FORMULA FOR Ap: 


From the tables of section 7.1 the numerical values of 
log A, according to 6.1.(65) are known for each value of P . 
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They prove to be almost constant for fixed values of p, 
wherefore the means of log A are taken and used. Consequently 


~ A, is considered to be independent of @ and only varying with 


The deviation for p = 100 is greatest for Q’, where it 
i. about 2.07%, and since no absolute accuracy is required the 
function Ay is sought in a form as simple as possible. 


As in 10 = 2.30259 we put 
(69) A_ = 1+ 2.30259 (log A_) (c + alog® +4e,), 
p P (p=100) P 
where c, a and b are constants. By putting & e," 0 this 
expression gives approximately 
Considering that A. as a maximum deviates only about 


st from unity we get, taking natural logarithms on both sides 
° (69) 


in A, = 2.30259 (log (c + a log® (B) +e) or 
log A 
= [c+ alog (P| + en 
(p=100) 


= c+ alog 

The left side of (70) are the relative values of the 
log A_’s measured by the value for p = 100; it will consequently 
be +120 for p = 100 and about -1.0 for p= 500. We determine 
the constants c, a and b by the requirement that e must have the 


correct values (i.e. eo" e,) for p = 20, 100 and 500. (Fig. 6). 
This gives the following formulae: 


(71) a = 1.023417 - 2.046836 (1 - 
(72) log b = 
(73) c = m- a log’ b 

where 


20) 


(74b) n = 4,00000 a - 0, 715338 (e, - 


and 


(74a) (° = 3.30103 a - 1.430677 (l-e 


= 
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(75a) m= + 1,30103 n - 1.69268 a 
(75b) M = €)99 + 2.00000 n - 4.00000 a 
(75c) M = + 2.69897 n 7.28444 
The formulae (74a) and (74b) must give the same value 


of n, and the three formulae (75) must give the same value of m, 
whereupon b and c are easily found by (72) and (73). 


Finally it follows from (69) with 
Ay = Cy + log” where 
(76b c, = 1 + 2.30259 c (log A 
P (p=100) 


(76c) a, = 2.30259 a (log A) 


0 (p=100) 


When the values of e° by 6.3.(70a) are deducted from the 
correct values of ep by 6.3.(70) we find & ep and consequently the 
relation between the errors from using the formulae (76a)b)c) and 
the errors produced for p = 100 by using the corresponding simple 
power formula. 

6.4. DERIVATION OF EXPRESSIONS FOR ay) b, AND 4, IN THE 
FORMULA 


1 - a, log 6) 
1 - a, log’ (F-) 
It will later be useful to have a formula like that above 
valid for values of by and by, which do not differ much. 


The values of a,, b, and a are to be expressed by 


the known quantities a), by, ay and bo. We find: 


= 1+4, - a, log’ 


a, log b, - a, log Do 


(79) log b, = - 
3 ay ao 


(80) 23° 


2 
(log b, - log bo) 
"1782, 1 - a+ log* (PY 
2 
) 
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In the practical application of (77), when b, and by, 
differ only slightly, 4, will often be insignificant, and the 
variation with p of the denominator [. - a log’ By. which is 
near 1, can be neglected, so that 


(80a) 4, (log b, - log by)’. 


6.5. DERIVATION OF FORMULAE FOR A@ AND ABSCISSAE 
§ OF THE FORM 


A@ = a+ bcot@ + + . 


The coefficients A@ for our different dimensionless 
quantities will be known for the four different values off . 


The same is true of the abscissae bj which are known 
for fixed values of n and the four values of @ . 


All these quantities can conveniently be expressed in 
formulae of the form (81), i.e. series with four members in 
increasing powers of the argument m = cot@. Since the inter- 
vals in @ are equal,’ the intervals in m will be unequal. To 
calculate the coefficients a, b, c and d so that the correct values 
of the function are obtained for @ = 15°, 20°, 25° and 30° , 
Newton’s Method of Interpolation is used. 


From the function 
(81a) g(m) = a + bm + cm® + dm® 
the values are known for the arguments Mo» My, My and m3 
(equal to cot 15°, cot 20°, cot 25° and cot 30°). 


We define: 
g(m) - g(mpo) 
(82) g,(m) = 
g,(m) - g,(m,) 
(83) g.(m) ——— = 
&_(m) - go(mo) 
(84) ga(m) = 


m- 


and calculate the special values g,(m,), &o(Mo) and g,(m,); this 
is done in the usual way in a table. 
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For the coefficients we then find the following expressions: 
(85) a= g(mo) + ma- g,(m,) + ml + - my 
(86) be @,(m,) - (mot mM + (mp my M>): g4(m,) 
(87) c¢ = go(my) - (mp + m, + 
(88) d= 


7. NUMERICAL DETERMINATION OF A@, Be AND 
A, FOR Q’, F’, v,,’. R’, (B+6o), AND v,’. 
* 
In the following sections the coefficients Ag, Boe and 
Ay for each of the dimensionless quantities Q’, F’, Vou e Tha 
(BB), I’ and Ya are denoted by the first subscripts gq, f, v, r, 


b, i and ti e.g. b » 
and respectively, for Q’ e.g Age ay and 


7 1. NUMERICAL CALCULATION OF Ag, Be AND log A,. 


m 


For the quantities Q’, F’, R’, (6+6,)) and I’ these 
calculations are made for the four values of + and for 
Q’ in Tables 7. 1/4a/b/c/d/ 


7.1/5a/b/c/d/ 
R’ - - 7.1/6a/b/c/d/ 
(B+6,) 7.1/%a/b/c/d/ 
I’ - 7.1/8a/b/c/d/. 
For Vin and ‘. the corresponding values of Ag, Be 


and log Ap can easily be found by combination of the above results. 


3 
| 
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af 


Q’ = 1. 244754 1n®(11 p) + 2.081904 + In (11 p) - 1.040952 


Table 7. 1/4a/ 


= 21.8184 Bog = 0.26247 log = 1.33882 
Q’ log logp Age* log A 
af Bo log P 
q? 
10 36.2472 1.55926 1.00000 0.26247 1.60129 -0.04201 
20 46,3995 1.66652 1.30103 0.34148 1.68030 -0, 01378 
$0 52.8926 1.72339 1.47712 0.38770 1.72652 -0, 00313 
50 61.6557 1.78995 1.69897 0.44593 1.78475 +0, 00520 
100 74,5852 1.87265 2.00000 0.52494 1.86376 +0.00889 
200 88.7109 1.94797 2.30103 0.60395 1.94277 +0. 00520 
300 97.5282 1.98913 2.47712 0.65017 1.98899 +0.00014 
500 109.2195 2.03830 2.69897 0.70840 2.04722 -0. 00892 
1000 126.1224 2.10080 3.00000 0.78741 2.12623 -0.02543 
= 21. 8184+ 26247 2 07% for 
24.6 500. 
ip = 20°) Table 7.1/4b/ 
= 0.974020 1n?(11 p) + 1.634774 In (11 p) - 0.817387 
Age 17.094 Bag 0.262349 log 1.23283 
B__-logp 
qe 
10 28,3873 1.45313 1.00000 0.26235 1.49518  -0.04205 
20 36.3354 1,56033 1.30103 0.34132 1.57415 -0, 01382 
30. 41.4186 1.61720 1.47712 0.38752 1.62035 -0, 00315 
50 48.2786 1.68376 1.69897 0.44572 1.67855 +0. 00521 
100 58.3999 1.76641 2.00000 0.52470 1.75753 +0, 00888 
200 69.4572 1.84171 2.30103 0.60367 1.83650 +0,00521 
300 76.3590 1.88286 2.47712 0.64987 1.88270  +0.00016 
500 85.5104 1.93202 2.69897 0.70807 1.94090 _-0, 00888 
For Ag, = 1: Q’ = 17.094 + + 2.07% for 


24.6 * p < 500. 
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Q’ = 0.805268 In*(11 p) + 1.355772 In(11p) - 0.677886 


Table 7.1/4c/ 


14,143 B__.= 0, 26230 log 1.15055 
a? qe “ae 
log 
@ log logp log A 
10 23.4869 1.37082 1.00000 0.26230 1.41285  -0,04203 
20 30.0609 1.47800 1.30103 0.34126 1.49181 -0,01381 
30 34,2651 1.53485 1.47712 0.38745 1.53800 -0.00315 
50 39.9387 1.60140 1.69897 0.44564 1.59619 +0.00521 
100 48.3094 1.68403 2.00000 0.52460 1.67515  +0,00888 
200 57.4539 1.75932 2.30103 0.60356 1.75411 +0.00521 
300 63.1617 1.80045 2.47712 0.64975 1.80030 +0.00015 
500 70.7297 1.84960 2.69897 0.70794 1.85849 -0.00889 
For 1: Q = 14,143 + 2.07% for 
24.6 £ p< 500. 
= Table 7.1/4d/ 
Q’ = 0.688370 In*(11 p) + 1.162312 In(11 p) - 0.581156 
12.0 = 0. log = 1.08276 
12.099 qe 0262249 log Ag, 
log 
Pp Q log Q logp Bag Bag log p log 
10 20,0915 1.30301 1.00000 0.26225 1.34501 -0,04200 
20 25.7135 1.41016 1.30103 0.34119 1.42395  -0,01379 
30 29.3087 1.46700 1.47712 0.38737 1.47013 -0,00313 
50 34.1605 1.53353 1.69897 0.44555 1.52831 +0,00522 
100 41.3183 1.61614 2.00000 0.52450 1.60726 +0,00888 
200 49.1376 1.69142 2.30103 0.60344 1.68620 +0.00522 
300 54.0182 1.73254 2.47712 0.64962 1.73238  +0.00016 
500 60.4893 1.78168 2.69897 0.70780 1.79056 -0.00888 
For Aj, =1: Q = 12.099- 792949 4 2.07% for 
24.6 500. 
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Table 7.1/5a/ 
F’= 1,244754 ln (11 p) + 2. 704281 


= 6.8585 = 0.109208 log = 0.83623 
Pp F’ logF’ logp Brg. log Arp + log Ay 
Brp- log p P 
3 10 8.55522 0.93223 1.00000 0.10921 0.94544 -0.01321 
20 9.41802 0.97396 1.30103 0.14208 0.97831  -0.00435 
, 30 9.92272 0.99663 1.47712 0.16131 0.99754 -0.00091 
50 10.55858 1.02360 1.69897 0.18554 1.02177 +0.00183 
100 11.42138 1.05772 2.00000 0.21842 1.05465 +0.00307 
200 12.28418 1.08935 2.30103 0.25129 1.08752 +0.00183 
300 12.78888 1.10683 2.47712 0.27052 1.10675  +0.00008 
500 13.42474 1.12791 2.69897 0.29475 1.13098  -0,00307 
1000 14,28754 1.15496 3.00000 0.32762 1.16385  -0.00889 
For A, = 1; F’= 6.8585 . 109208 9 
for 23.7 3 500. 
F’ = 0.974020 In(11 p) + 2.121784 
Agg = 5.3728 = 0.10911 log Ayg, = 0.73020 
pF log  logp log peAret og Ag 
- log p 
fp 
10 6.70015 0.82608 1.00000 0.10911 0.83931 -0.01323 
20 7.37529 0.86778 1.30103 0.14195 0.87215  -0.00437 
30 7.77021 0.89043 1.47712 0.16117 0.89137 --0.00094 
50 8.26777 0.91739 1.69897 0.18537 0.91557  +0.00182 
100 8.94291 0.95148 2.00000 0.21822 0.94842 +0.00306 
200 9.61806 0.98308 2.30103 0.25106 0.98126 +0.00182 
300 10.01298 1.00056 2.47712 0.27027 1.00047 +0.00009 
500 10. 51054 1.02162 2.69897 0.29448 1.02468  -0.00306 
For A, = 1: F’= 5.3728. 9. 707% 
'P for 238£p < 500 
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F’ = 0.805268 In(11 p) + 1. 758406 


Table 7.1/5c/ 


57 


4.4456 Brg, = 0.10906 log Arg =0. 64793 
Pp F’ log F’ logp Br, - log ps Argt log As 
? B log p P 
{p° 
10 5.54355 0.74378 1.00000 0.10906 0.75699  -0.01321 
20 6.10172 0.78545 1.30103 0.14189 0.78982 -0.00437 
30 6.42823 0.80809 1.47712 0.16109 0.80902 -0.00093 
50 6.83958 0.83503 1.69897 0.18529 0.83322 +0.00181 
100 7.39775 0,86910 2.00000 0.21812 0.86605 +0.00305 
200 7.95593 0.90069 2.30103 0.25095 0.89888  +0.00181 
300 8.28243 0.91816 2.47712 0.27015 0.91808  +0.00008 
500 8.69379 0.93921 2.69897 0.29435 0.94228  -0.00307 
For A. =1: F’= 4.4456. 19906 4 9 70% 
for 238% p < 500. 
Table 7.1/5d/ 
F’ = 0.688370 In (11 p) + 1.506497 
Arg = 3.8034 Brg = 0.10899 log Arg = 0.58017 
P F’ log F’ log p Bro -log p lofAre + log Arp 
Bro - 108 P 
10 4.74217 0.67598 1.00000 0.10899 0.68916 -0.01318 
20 5.21931 0.71761 1.30103 0.14180 0.72197 -0.00436 
30 5.49842 0.74023 1.47712 0.16099 0.74116 -0.00093 
50 5.85006 0.76716 1.69897 0.18517 0.76534 +0.00182 
100 6.32720 0.80121 2.00000 0.21798 0.79815  +0.00306 
200 6.80434 0.83278 2.30103 0.25079 0.83096 +0.00182 
300 7.08345 0.85025 2.47712 0.26999 0.85016 +0.00009 
500 7.43509 0.87129 2.69897 0.29417 0.87434 -0.00305 
For A, = 1: F’ = 3.8034. 4 0. 107% 
fp for 23.8 $p< 500. 
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Table 7.1/6a/ 


R’ 1, 244754 In(11 + 2.704281 
. n p + . 


0.73272 0.020945 log A, = 9.86494 
log + 
R’ log R’ logp Bip log P log Arp 
rp .log p 
10 0.760165 9.88091 1.00000 +0.02095 9.88589 -0.00498 
20 0.777242 9.89056, 1.30103 +0.02725 9.89219 -0,00163 
30 0.786149 9.89550 1.47712 +0.03094 9.89588 -0.00038 
50 0.796406 9.90113. 1.69897 +0.03558 9.90052 +0.00061 
100 0.808846 9.90786 2.00000 +0.04189 9.90683 +0.00103 
200 0.819853 9.91374 2.30103 +0.04820 9.91314 +0.00060 
300 0.825724 9.91683 2.47712 +0.05188 9.91682 +0.00001 
500 0.832597 9.92043 2.69897 +0.05653 9.92147 -0.00103 
1000 0.841099 9.92485 3.00000 +0.06284 9.92778 -0,00293 
For A, = 1: = 0.73272. "#0845, 0,237 % for 24. 8Sp< 500. 
220°, Table 7,1/6b/ 


R’ - 0.974020 In (11 p) + 2.121784 
n Pp + 


0.73172 0.02103 log 9. 86434 

p log R’  logp log pe “re* log 

rp og p 

10 0.759273 9.88040 1.00000 0.02103 9.88537 - 0.00497 
20 0.776382 9.89008 1.30103 0.02736 9.89170 -0.00162 
30. 0.785307 9.89504 1.47712 0.03106 9.89540 -0.00036 

50 0.795586 9.90068 1.69897 0.03573 9.90007 +0.00061 
100 0.808056 9.90744 2.00000 0.04206 9.90640 +0.00104 
200 0.819092 9.91334 2.30103 0.04839 9.91273 +0.00061 
300 0.824979 9.91644 2.47712 0.05209 9.91643 +0.00001 
500 0.831871 9.92006 2.69897 0.05675 9.92109 - 0.00103 
For A, = 1: R’= 0.73172 . 92103 240 % for2h6! p < 800. 
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R’ 0.805268 . ln (11 + 1, 758406 


rp 


0.72968 


Table 7.1/6c/ 


r@ 


0.02118 


log 


= 9.86313 


59 


P 


R’ 


log R’ 


log p 


B rp 8 


P 


log A 


+ 
re log Arp 


10 
20 
30 
50 
100 
200 
300 
500 


0.757292 
0. 774488 
0. 783462 
0. 793800 
0.806345 
0.817451 
0.823376 
0.830315 


9,87927 
9.88901 
9.89402 
9.89971 
9.90652 
9.91246 
9,91560 
9.91925 


1, 00000 
1, 30103 
1.47712 
1.69897 
2.00000 
2.30103 
2.47712 
2.69897 


0.02118 
0.02755 
0.03128 
0.03598 
0.04235 
0.04873 
0.05246 
0.05716 


9. 88431 
9.89068 
9.89441 
9.89911 
9,90548 
9.91186 
9.91559 
9.92929 


- 0,00504 
- 0,00167 
- 0.00039 
+0, 00060 
+0,00104 
+0,00060 
+0,00001 
- 0.00104 


For Arp = R’= 0.72968 . p 


0.02118 + 9 240 % for24.9$ pz 500. 


A 
rp 


0. 72638 


Table 7.1/6d/ 


R’ = 0.688370 . In (11 p) + 1.506497 
. n Pp + 


Pre 


= 0,021460 


log A. 


= 9.86117 


R’ 


log R’ 


log p 


Pa 


log A 


+ 


log A 


rp 


10 
20 
30 
50 

100 

200 

300 

500 


0. 754333 
0.771664 
0.780713 
0.791141 
0.803801 
0.815014 
0.820998 
0, 828008 


9.87756 
9.88743 
9.89249 
9,89825 
9.90515 
9.91117 
9.91434 
9.91803 


1.00000 
1, 30103 
1,47712 
1,69897 
2.00000 
2.30103 
2.47712 
2.69897 


0.02146 9.88263 
0.02792 9.88909 
0.03170 9,89287 


0.03646 
0.04292 
0.04938 
0.05316 
0.05792 


9.89763 
9.90409 
9.91055 
9.91433 
9.91909 


- 0.00507 
- 0.00166 
- 0,00038 
+0,00062 
+0.00106 
+0,00062 
+0.00001 
- 0.00106 


For A. 


zis 
P 


R’ = 0.72638. p 


0.021460 


+ 0.244 % for24.7£ p < 500. 
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Table 7,1/7a/ 


B+ Bo = 1, 244754 In (11 p) + 5.318022 


= 9.2804 


Bho 


= 0.088795 


log A 


bp 


= 0.96757 


log 


+ 


log A 
- lo b 


50 
100 
200 
300 
500 

1000 


11. 16896 
12,03177 
12. 53646 


13.17232 


14,03512 
14, 89792 
15. 40262 
16, 03848 
16,90128 


1.04801 
1.08033 
1,09818 
1.11967 
1.14722 
1.17313 
1.18759 
1, 20516 
1, 22792 


1, 00000 
1, 30103 
1,47712 
1.69897 
2.00000 
2.30103 
2.47712 
2.69897 
3.00000 


0.08880 
0.11552 
0.13116 
0.15086 
0.17759 
0, 20432 
0, 21996 
0, 23966 
0. 26639 


1, 65637 
1, 08309 
1,09873 
1.11843 
1,14516 
1,17189 
1,18753 
1, 20723 
1, 23396 


- 0.00836 
- 0.00276 
- 0.00055 
+0,00124 
+0,00206 
+0.00124 
+0, 00006 
- 0.00207 
- 0.00604 


For A,,)= 13; B+ Bo = 9.2804 . p 


0.088795 


+ 0.475 % for 23. 2$p< 500, 


3 


= 7.2375 


Table 7.1/7b/ 


B+ By = 0.974020 In (11 p) + 4.135563 
B, = 0.08901 
be 


log = 0.85959 


B+ Bo 


log(B+By) log p 


bp 


+ 


30 


8.713925 


9. 389066 
9.783993 


50 10,281551 
100 10.956693 
200 11.631835 
300 12.026761 
500 12.524320 


0.94022 1.00000 


0.97262 
0.99052 
1.01206 
1, 03968 
1.06565 
1.08015 
1.09775 


1, 30103 
1.47712 
1.69897 
2.00000 
2.30103 
2.47712 
2.69897 


0.08901 
0.11581 
0.13148 
0.15123 
0.17802 
0.20482 
0, 22049 
0.24024 


0.94860 
0.97540 
0.99107 
1.01082 
1.03761 
1.06441 
1, 08008 
1.09983 


- 0.00838 
- 
- 0,00055 
+0.00124 
+0,00207 
+0.00124 
+0.00007 
-0,00208 


For A 


pp 


B+B, = 7.2375. p 


0.08901 


+ 0.477% for23.2£p < 500 


i 60 
? = 
: 
10 
20 
30 
: 
4 log A | 
log p p 
10 
|_| 


B+ By = 0.805268 ln (11 p) + 3.400201 


Table 7,1/7c/ 
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Ang = 5.9660 Bug = 0.08919 log Ano = 0.77568 
log Angt 
p B+ Bo log(B+ Bo) log p By glog p log 
10 7,185347 0.85645 1.00000 0.08919 0.86487 -0.00842 
20 =7.743519 0.88894 1.30103 0,11604 0.89172 -0.00278 
30 8.070023 0.90687 1.47712 0.13175 0.90743 -0. 00056 
50 8.481378 0.92847 1.69897 0.15154 0.92722 +0,00125 
100 9.039549 0.95615 2.00000 0.17839 0.95407 +0.00208 
200 9.597721 0.98217 2.30103 0. 20524 0.98092 +0,00125 
300 9.924225 0.99670 2.47712 0.22094 0.99662 +0,00008 
500 10.335580 1.01433 2.69897 0.24073 1.01641 -0.00208 
For Ay, = 1: = 5.9660 + 0, 48% for23.2Sp 500, 


ip = 30°, Table 7.1/7d/ 


B+ By = 0.688370 In (11 p) + 2.891576 


Ang = 5-0864 = 0.08936 log 0.70641 
log A et 

P log(6+B)) log p Big: 10g P Byglog log Any 
10 6.12774 0.78730 1.00000 0.089360 0.79577 -0.00847 
20 6.60488 0.81986 1.30103 0.116260 0.82267 -0,00281 
30 6.88399 0.83784 1.47712 0.131995 0.83841 -0,00057 
50 7.23562 0.85948 1.69897 0.151820 0.85823 ¥+0.00125 
100 7.71277 0.88721 2.00000 0.178720 0.88513 +0,00208 
200 8.18992 0.91328 2.30103 0.205620 0.91203 +0.00125 
300 8.46902 0.92783 2.47712 0.221355 0.92777 +0,00006 
500 8.82067 0.94551 2.69897 0.241180 0.94759 -0.00208 
For 1: 6+ = 5.0864. p? 98936 48% for23.3$ p z 500, 


Me 
& 
| 


Avo 


= 6.9199 


Table 7.1/8a/ 
I’ = 1, 244754 (11 p) + 2.771323 


Bip 


= 0.108561 


log A, 


= 0.84010 


? 


log I’ 


log p Bio: log p B. 
i 


log A 


ie * 
Q' log p 


log Aip 


100 
200 
300 
500 
1 000 


8.62226 
9.48507 
9.98976 
10.62562 
11, 48842 
12,35122 
12, 85592 
13.49178 
14,35458 


0.93562 1.00000 


0.97704 
0.99957 
1,02635 
1, 06026 
1.09171 
1.10910 
1,13007 
1.15699 


1.30103 
1.47712 
1,69897 
2.00000 
2.30103 
2.47712 
2.69897 
3.00000 


0.10856 
0.14124 
0.16036 
0.18444 
0.21712 
0, 24980 
0.26892 
0, 29300 
0.32568 


0.94866 
0.98134 
1, 00046 
1,02454 
1.05722 
1, 08990 
1,10902 
1.13310 
1, 16578 


- 0.01304 
- 0, 00430 
- 0, 00089 
+0, 00181 
+0. 00304 
+0, 00181 
+0, 00008 
- 0.00303 
- 0, 00879 


For A; 


=1; I’ = 6.9199. p 


0.108561 


+ 0.702% for 23.69% p< Soo. 


le - 20°, | 


I’ = 0,974020 In (11 p) + 2, 207429 


Table 7,1/8b/ 


Aig = 5.4502 Bye = 0,108096 log = 0.73642 
p I log I log p Bio: log p Bip" log p log Aiy 
10 6.78579 0.83160 1.00000 0.10810 0.84452 -0,01292 
20 7.46093 0.87280 1.30103 0.14064 0.87706 -0,00426 
30 7.85586 0.89519 1.47712 0.15967 0.89609 -0,00090 
50 8.35342 0,92186 1.69897 0.18365 0.92007 +0.00179 
100 9.02856 0.95562 2.00000 0.21619 0.95261 +0.00301 
200 9.70370 0.98694 2.30103 0.24873 0.98515 +0,00179 
300 10,09863 1.00426 2.47712 0.26777 1.00419 +0.00007 
500 10,59619 1.02515 2.69897 0.29175 1.02817 -0.00302 
For A, 1:1" = 5.4502, 9. 695% for 23, 72 pz 500, 


: 
c 
62 
: 
i 


I’ = 0,805268 In (11 p) + 1.862054 


A; = 4.5399 


Table 7.1/8c/ 


= 0.107547 


log A; 


= 0.65705 
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P 


log I’ 


log A. 


+ 
log p Bip P p log A; 


10 
20 


5.64720 
6, 20537 
6.53188 
6.94323 
7.50140 
8.05957 
8. 38608 
8.79743 


0. 75183 
0.79277 
0.81504 
0, 84156 
0.87512 
0.90631 
0.92355 
0.94435 


1, 00000 
1, 30103 
1.47712 
1.69897 
2.00000 
2.30103 
2.47712 
2.69897 


0.10755 
0.13992 
0.15886 
0.18272 
0.21509 
0.24747 
0. 26641 
0, 29027 


0. 76460 
0, 79697 
0.81591 
0.83977 
0.87214 
0.90452 
0.92346 
0.94732 


- 0,01277 
- 0,00420 
- 0,00087 
+0.00179 
+0.00298 
+0.00179 
+0. 00009 
- 0.00297 


For A; 


: I’ = 4.5399 p 


0.107547 + 0.688% for 23.70 f p< 


500, 


" 


0.688370 


= 3.9143 


Table 7, 1/8d/ 


In (11 p) + 1.627882 


Bip 


= 0.106933 


log A; 


= 0.59266 


log I’ 


log A. + 


log p Bip log p Big: p log Aip 
i 


50 
100 
200 
300 
500 


4.86355 
5. 34070 
5.61980 
5.97144 
6.44859 
6.92573 
7.20484 
7.55648 


0.68695 
0.72760 
0.74972 
0. 77608 
0. 80947 
0. 84046 
0.85762 
0.87832 


1,00000 
1, 30103 
1.47712 
1.69897 
2.00000 
2. 30103 
2.47712 
2.69897 


0.10693 
0.13912 
0.15795 
0.18168 
0, 21387 
0, 24606 
0. 26489 
0, 28861 


0.69959 
0.73178 
0.75061 
0.77434 
0.80653 
0.83872 
0.85755 
0.88127 


- 0.01264 
- 0.00418 
- 0,00089 
+0.00174 
+0.00294 
+0.00174 
+0.00007 
- 0,00295 


for A. 
i 


P 


3.9143. p 


0.106933 


+ 0.679% for 23.76 * p< 500. 


= 
30 
50 
100 
200 
300 
500 
10 
20 
| 30 
| 
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7.2. NUMERICAL CALCULATIONS OF FORMULAE FOR Ay AS 
FUNCTION OF p. 


These calculations are compiled in Tables 
7.2./9/10/11/12/ for 
Aap Ary Abp Aip » respectively. 
The values of log Ay are found in the last columns of the 
tables of section 7.1. It is seen that the values of e.g. A p are 
practically equal for the different values of P, and the same 


applies to the other A, s. The values in the tables of section 7.2. 


are taken as the mean values and denoted by (log A)) 


The corresponding values of ep from formula 6.3.(70) are 
next found, and the values C09 ©100 and Cs are introduced in 


formulae 6.3.(71 through 75), which give the values of a, b and c 
indicated in the following tables. These values introduced in 


formulae 6.3.(76a)b)c) thus give the definite formulae for the A,'s. 


The errors of these formulae are found by calculation of 
the values of ep according to 6.3.(70a) 


eo = ¢ + a log*(P). 


P 


The differences ey - ep 2 ae, are proportional to the 


errors; the absolute values of the aaae are found by multi- 
plication by ae, of the errors for p = 100 stated in the tables of 


section 7.1. 


F 
| 
| 
4 
§ 
: j 
a 


| 
3 
3 


A... Table 7.2/9/ 
p (log 6.3(70a) 4ep 
10 -0.04202 -4,.7320 -4.0575 -0.6745 
20 -0.01380 -1.5541 0.0000 
30 -0.00314 -0.3536 -0.4815 +0.1279 
50 +0,00521 +0.5867 +0.4585 +0.1272 
100 [+0.00888] [+1.0000] +1.0000 0.0000 
200 +0.00521 +0.5867 +0.6967 -0.1100 
300 +0.00015 +0,0169 +0.1225 -0.1056 
500 -0.00889 -1,0011 0.0000 
1000 -0,.02543 -2,.8637 - 3.2659 +0.4022 


By 6.3 (71 through 75) : 


a =-4,6619 , b= 110.262, c = 1.0084 , log b = 2.04243 
By 6.3(76a)b)c): for ~ 17.5 ~ 622 
Aap = 1.020619 - 0.095322 0.265 % 
Ary: Table 7. 2/10/ 
by 
p loge), 6.3(70a) 
10 -0,01321 -4.3170 -3.8322 -0.4848 
20 -0.00436 -1.4247 -0.0001 
30 -0.00093 -0.3039 -0.3968  +0.0929 
50 +0,00182 +0.5948 +0.4980 +0.0965 
100 [+0,00306] [+1.0000) +1.0002 +0.0002 
200 +0.00182 +0.5948 +0.6813 -0.0865 
300 +0.00009 +0.02941 +0.1143 -0.0849 
500 -0,00306 -0.9997 -0,0003 
1000 -0.00889 -2.9053 -3.2243 +0.3190 
By 6.3 (71 through 75) 
= -4,52842 , b = 108.035, c = 1.0053, logb = 2.03354 
By 6.3(76a)b)c): for 17.4 £ p< 616 
Ag, = 1.007083 - 0.031907 - log” 0.068 % 


65 
| | 
be 


App: Table 7.2/11/ 
Cp = ep by 
10 -0.00840 -4,0508 -3,6917 -0,3663 
20 -0,00278 -1,3431 +0,0001 
30 -0,00056 -0,.2705 -0.3432 +0,0727 
50 +0,00125 +0,6039 +0.5238 +0.0801 
100 [+0.00207] [+1.0000] +0.9998 +0.0002 
200 +0.00125 +0.6039 +0,6694 -0.0655 
300 +0.00007 +0,0338 +0,1023 -0.0685 
500 -0,00208 -1,0050 +0,0002 
1000 -0.00604 -2.9179 -3.2079  +0,2900 
By 6.3 (71 through 75) : 
a=-4.44962 , b= 106.46 , c = +1.0031 , log b = +2,02718 
By 6.3 (76a)b)c): for17.2 p< 605 
= 1-0047812 - 0.0212087 - log” 0.0382 % 


Table 7. 2/12/ 


10 -0.01284 -4.2943 -3,8208 -0,4735 
20 -0.00424 -1,4182 +0,0001 
30 -0.00089 -0.2977 -0,3927 +0.0950 
50 +0,00179 +0.5987 +0.5001 +0,0986 
100 [+0,00299] |+1,0000) +0.9999 +0,0001 
200 +0.00178 +0.5953 +0.6801 -0,0848 
300 +0.00008 +0.0268 +0.1132 -0,0864 
500 -0.00299 -1,0002 +0,0002 
1000 -0.00879 -2.9398 -3.2229 +0. 2829 


By 6.3 (71 through 75) : 


a = -4.5216 , b= 107.914, c = +1.0048 , logb = +2.03307 
for 17.3 £ p < 636 


+ 
ip 1-0069177 - 0.031130 log* 0.0679 % 


By 6.3 (76a)b)c) : 


A 


| 
| 
| | 
| 
ip 
: 
‘ 


SEC. 8.1. 67 


8. GENERAL REVIEW OF FORMULAE FOR Q’, 
F’, v.,’, R’, B’, I’ AND v,’ AND OF 
FORMULAE FOR THE ABSCISSAE § AS 


FUNCTIONS OF @ FOR FIXED VALUES OF 
THE ORDINATES 1 - 


8.1. THE DEFINITIVE FORMULAE FOR THE DIMENSIONLESS 


QUANTITIES Q’, F’, Yn’ B’, I’ AND 


These are given in the following tables: 
Q’ in Table 8,1./13/, Formulae (89) and(89a)b)c)d)e) 
F’ mo . (90) - (90a)b)c)d)e) 
Vin’ - - 8.1/15/, (91) - (91a)b)c)d)e) 
- - 8.1./17/, (93) - (93a)b)c)d)e) 
I’ o - (94) - (94a)b)c)d)e) 
- _63./10/, (95) - (95a)b)c)d)e) 


The coefficients Age etc. are in general functions of @ ; 


they are given as functions of m = cot@; in some cases the 
variation with@ is so slight that a mean value near the correct 


value for @ = 20° is used. 
The exponents B@ in most cases vary so little with @ 


that a mean value can be used. The factors A etc. are 


generally functions of p = — as indicated. For R’, B’ and 
Val these factors can be put equal to unity, because the simple 
power formulae have already very small errors. 

In the tables are given the directly calculated values for 
Ag and Be for @ = 15°, 20°, 25° and 30°. 

The tables also state the ranges of validity of the simple 


power formulae and of the complete formulae, including the 
factors A, with the corresponding maximum errors as well. 


j 
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Table 8.1/13/ 


B 
| Age 
15° | 21.8184| 0.26247 
20° | 17.094 | 0.26235 
25° | 14.143 | 0.26230 
Q for half of total section. 30 12.099 0. 26225 


(89c) Ago = 3.1815. Age 


8 = 0.26234 
(89d) Bue mean 262 


(89e) = 1.020619 - 0.095322-1og* 


For A, = 1 and 25 £ p < 500; 
B 


(89a) Q’ = P + 2.07 % 


Complete formula and 17.5 £ p 5 622: 


(89) 
B 
(89b) Q =A .p +t 0.265% 
a? ap 
} 


Table 8.1/14/ 


F 
B 
max” 
15° |6.8585 |0.10921 
p 7. A 20° | 5.3728 |0.10911 
fp 
25° | 4.4456 |0.10906 
F for half of total section] 30° | 3.8034 |0. 10899 


fp 


(90c)  Agg= 1.03153 + 1.63688-m - 0.02187-m? + 0.000437-m® 
(90d) Brg mean = 0.1091 
(90e) Ap, = 1.007083 - 0.031907 . log (rob-wss) 
For Ay, = 1 and 25 £p < 500; 
B 
(90a) F’ =A, p @t 0.71% 


B 
(90b) F’ = Age. p 


Complete formula and 17.4 $p<z 616. 


+ 


6 
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Table 8.1/15/ 


V, 
g 15> | 3.1815 | 0.15324 
A 20° | 3.1816 | 0.15324 
25° | 3.1813 | 0.15324 
30° | 3.1811 | 0.15326 


Q9lc) A = 3.1815 
( 91c) ve mean 


91d) B => 0.15324 
(91d) vp meen 


(91e) AY, = 1.0135645 - 0, 0634373+ log” 


For A, = 1 and 25 p < 500; 


B 
ve + 
(91a) Ave’ p 1.37 % 


Complete formula and 17.5 pz 620 
B 
(91b) =A. p Ay 0.195 % 


‘2 
q 
# 
- 
: 
4 


Table 8.1/16/ 


R 
(02) A B 
re re 
15° | 0.73272 | 0.02095 
20° | 0.73172 | 0.02103 
25° | 0.72968 | 0.02118 
30° | 0.72638 | 0.02146 
(92c) A, mean > 0.732 for15 20° 
(92d) BL, mean = 0.0210 - 182492 20° 
(92e) — = 1 as simple power-formula sufficient 
For = and 25 p< 500 
B 
(92a,b) R’ =A. . + 9.24% 


re? 


71 


q 
4 
q 
4 
i 


72 


Table 8.1/17/ 


4 s B s 
(93) BY = 5 (B + Bo) ? | Ane Bie 
Bug 15° | 9.2804 | 0.08880 
20° 7.2375 | 0.08901 
25° | 5.9660 | 0.08919 
B for half of total section. 30° 5.0864 | 0.08936 


(93c) Ang 1,29443 + 2.23919'm - 0.03069-m? + 0. 00109-m" 
(93d) Bue mean = 0.0890 


(93e) Ay, = 1.0047812 - 0.0212087-log” 


For Ay, = 1 and 25 £p < 500 


B 
(93a) BY = p P+ 0.48% 


la 


Complete formula and 17.2 p < 605 


93b) = Ang - P Ay, = 0.038 % 


( 
ES 
° 
| 
| 
4 


wat. 


Table 8.1/18/ 
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(94) I’ = _g-F.1 A, B, 
Tmax 2 
Ymax'—9 15° | 6.9199 | 0.108561 
20° | 5.4502 | 0.108096 
A 
ig ip 25° | 4.5399 | 0.107547 
30° | 3.9143 | 0.106933 


(94c) = 1.27988 + 1.51750-m + 0.00519-m* - 0. 00184-m* 


(94d) B, mean = 0,10833 for 15°S p 2 20° 


(94e) A, = 1.0069177 - 0.031130-log* 


‘ 
For Ain = } and 25 = p z 500 


B 
(94a) =A. p 1? t 9.70% 


1? 


Complete formula and 17.3 p < 636 


B. 
+ 


(94b) I’ = Aig ip 0.068 %. 


q 
+ 
4 
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Table 8,1/19/ 


Aye Bue 
15° | 1.1630 | - 0.01015 
B 
20° 1.1604 | - 0.00997 
*P | 25° | 1.1584 | - 0.00983 


30° 11.1566 | - 0.00970 


(95c) mean 


xp 


(95e) Aye = 1 as simple power-formula sufficient. 


1.160 for 15° 2 25° 


0.0100 for 15° £P 25° 


82 


(95d) B mean 


ForA._ = 1 d 25 £ pe 500 
or A, and 25 5 


B 
(95a)b) = Aye .p 


| 
(95) v? 
4 


SEC. 8.2. 75 


8.2. FORMULAE FOR THE ABSCISSAE $ AS FUNCTIONS OF 
m= cot@ . 


These formulae are found by the method indicated in 
Sec. 6.5. and are based on the values of § , which are found 
for fixed values of ) and for the four different values of @ 
as listed in Tables 3.3. /2a/b/c/d/. 


The numerical coefficients a, b, c and d in the formulae: 
(96) S$ m+tcm’ + dm’, 


(96a) where m = cot@. 
are given in Table 8.2./20/. 


A separate formula for s is consequently available for 
each of the fixed values of 1). 


For later use we have computed by means of these 
formulae the numerical values of @ for q@ = 14°, 16°, 17°, 


18° and 19°". These are, together with the previously found values 


for = 15°, 20°, 25° and 30° (Tables 3,3, (2a)b)c)d) ), compiled 
in the Table 8.2./21/. 


4 
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Table 8,2/20/ 


n § 1° m: 

0 &, = 8 0 0 0 

0.1 $5.1 +0.00723 +0,10041 +0,00049 0, 000065 
0.2 +0,02808 +0.20256 +0.00179 - 0.00024 

0.3 +9-06579 +0,30833 +0.00365 - 0.00051 

0.4 $o.4 +0.11622 +0.42315 +0.00455 - 0.00070 

0.5 $o.5 +9-18331 +0.54901 +0.00447 0.00082 

0.6 $4 +0.27246 +0.68819 +0.00360 -0,00090 

0.7 +0.38875 +0,84935 +0.00049 - 0.00078 

0.75 +0.46333 +0.93745 - 0.00083 - 0.00078 

0.8 $, = +9.54583 +1,04230 -0,00555 - 0.00039 
0.85 $5 ys" +0.64936 +1.15793 - 0.01000 - 0.00007 

0.9 +0. 78211 +1.29522 -0.01590 +0. 00039 

0.95 $9 +0.97200 +1.47634 0.02432 +0.00107 

0.98 gg= +1.15392 +1.63955 - 0.03239 +0. 00175 

1.0 = +1,46574 +1.90032 -0.04441 +0.00274 


aS 
‘ 
: 
: + 
. 
— 


Table 8,2/21/ 


p= 14° p= 15° p= 16° 179 18° 
0 0 0 0 0 0 

0.1 0.41364 0.38543 0.36060 0.33863 0.31901 
0.2 0.85471 0.79740 0.74698 0.70227 0.66235 
0.3 1,32825 1.24081 1.16282 1.09549 1.03443 
0.4 1,84141 1.72240 1.61757 1.52446 1.44123 
0.5 2.40426 2.25187 2.11750 1,99816 1.89142 
0.6 3.03248 2.84419 2.67807 2.53045 2.39835 
0.7 3.75287 3.52483 3.32366 3.14480 2.98468 
0.75 4.15957 3.90979 3.68943 3.49342 3.31790 
0.8 4.61183 4.33820 4.09672 3.88201 3.68976 
0.85 5.12818 4.82792 4.56295 4.32735 4.11634 
0.9 5.74634 5.41474 5.12223 4.86212 4.62914 
0.95 6.57108 6.19867 5.87020 5.57814 5.31653 
0.98 7.32166 6.91261 6.55199 6.23136 5.94415 
1.0 8.54989 8.08172 7.66902 7,30216 6.97354 


717 


To be continued 


. 
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Table 8.2/21/ 


s 


19° p= 20° p= 25° 30° 
0 0 0 

0.1 0.30138 0.28545 0.22415 0.18226 
0.2 0.62648 0.59404 0.46923 0.38394 
0.3 0.97955 0.92989 0.73877 0.60813 
0.4 1.36637 1.29864 1.03770 0.85917 
0.5 1.79535 1.70843 1.37312 1.14338 
0.6 2.27942 2.17176 1.75597 1.47056 
0.7 2.84046 2.70984 2.20476 1.85730 
0.75 3.15977 3.01649 2.46225 2.08050 
0.8 3.51653 3.35955 2.75167 2.33244 
0.85 3.92618 3.75381 3.08587 2.62459 
0.9 4.41914 4.22878 3.49046 2.97984 
0.95 5.08068 4.86682 4.03673 3.46169 
0.98 5.68520 5.45034 4.53826 3.90564 
1.0 6.67722 6.40843 5.36379 4.63820 


(continued) 


} 
| 
| 
: 
Met 
a 


inued) 


SEC, 9. 79 


9. DETERMINATION OF ISOVELS IN 
THE EQUILIBRIUM PROFILE. 


The principal assumptions made in the foregoing sections, 
which made it possible to determine a definite equilibrium profile 
for fixed values of 9 , € maw F , k and Q, are expressed by the 


formulae Ah 2.1(3), 2.2(9), 4.2.(19), 4.4.(45), giving the 
distribution of the shearing stresses, the relation between hydraulic 
lifting force and shearing stress, the logarithmic law of velocity 
distribution and the principle of minimum cross section. 


On the same assumptions it is possible to calculate the 
accurate form of the isovels and thus, by comparing these with the 
results from actual measurements, to get a certain check on the 
correctness of the assumptions. 


We therefore indicate the method of determination of 
isovels having velocities expressed as multiples of the mean 
velocity of the total cross section. 


The velocity at a point P of the cross section is denoted 
by ¥. . Here the subscripts indicate that P is situated at: the 


distance z from the bottom ina normal to this with the length 
& * Venew between the bottom and the water surface. By 4, 2,(18) 


and 4.2.(19) we find this velocity to be 


(97) Ve . [3.392 + 


Substituting 
(98) v 
sr and 
m 
(99) 
max 


y 
we get, with ai = p as usual: 


(100) or Ve - In [29. 725 p- 


Further we have from Table 8.1. /15/ 


Vv 
8.1.(91b)c)d) —™ 3.1815 15324. 4 
2.5\| max 


i 
4 
ia 


and 
8.1.(9le) = 1.013565 - 0.063437 log*( whence 


(01) or in| 20.725 p 
| 


.1815 p Avp 


This formula gives the velocity at the point P measured 
with the mean velocity Yn 2&8 unity as function of the dimension - 


less quantities & , total length of the normal, z’, distance from 


the bottom, and p = 


If we choose the normals a used in the numerical 


integration, these will vary with @ ; but if the normals are chosen 
as fixed multiples of , introduced in a drawing of the 


equilibrium profile in question, the relative velocities will depend 
only on p, i.e. the relative roughness of the bottom. 


From (101) we get: 


3.1815 15324 | Ave 
29.725: p 


For a fixed value of p a complete set of isovels, e.g. 
with relative velocities 
can be drawn by choosing suitable values of © . 


As an example the isovels are worked out for an equi- 


librium profile with @ = 30° and p = ex - 50. 
We get 

29.725 » p = 1486.25 and from 8.1. (91e) 
A 1.00588 


vp 
3.1815 pr 155% . 5.8281 
5.8281 
r) 
(10laa) 2° = 


For the values § we have chosen the values Cm for 


which V5m can be found in Table 4.6. /3d/. 


ae 80 SEC. 9. 
| 
| 
i | 
— 
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The results are shown in Fig. 7. 


It will be seen that the form of the isovels compares 
fairly well with the results of actual measurements, and it can 
consequently be concluded that the assumption made concerning 
the distribution of shearing stresses must be at least approx- 
imately correct. 


10. STUDY OF MODELTESTS CARRIED OUT 
IN VIENNA,1916. COMPARISON WITH THEORY. 


10.1. DESCRIPTION OF MODEL TESTS AND RESULTS. 


These model tests e. briefly mentioned in [2] 3 the 
complete report is found in | 3 

The sand used had a median grain diameter of about 

d 
1.4 mm and Ms = 1.8; it was a rather sharp quartz sand 
10% 

with 38% of voids, measured in loose filling, and with a specific 
weight of 1.64 kg/1. 


Three different profiles were shaped by three different 
constant discharges, and when the profile had become stationary, 
its size and shape was carefully measured by a photogrammetric 


method. The limiting tractive force for the sand, ‘max’ Was 


found by separate tests on a nearly horizontal bottom and turned 
out to be 0.075 kg/m*. 


For each equilibrium profile the discharge and the slope 
of the water surface (at uniform movement) were further measured. 


The observed results are given in the following 
Tables 10.1 /22/and 10.1 /23/. 


Table 10.1 /22 
Dimensions, Slope and Discharge, 


Width of Slope Area of Maximum | Discharge 
Test | water surface I cross section depth 2Q 
No 2B Too 2F y 
(cm) (cm?) tami (1/sec) 
1 100 1.35 446.2 5.6 13.8 
2 198 0.79 1444.0 9.0 53.5 
280 0.65 2475.0 10.5 87.0 


4 
a 
n 
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fable 10.1, /23/ 
Co-ordinates x and y. (See Fig. 3) 


Test | Depths Abscissae x 
No y em) 


y (1.6 9.8 3.8 4.6 6.29° 7.9 6.3°°8.9 8.2 


Test | Depths Abscissae x 
No y cm 
y 


2 y 8.6 8.8 8.9 9.0 9.0 
y 9.4 9.6 9.8 10.0 10.1 10.2 10.3 10.4 10.5 


10,2. COMPARISON OF THEORY WITH MODEL TESTS; 
COMMENTS. 


If we compare the results of these tests with the theory 
developed in the present paper, we should note that the area 
of the cross section necessary to carry a definite discharge 
without erosion varies only insignificantly with the angle of 
friction @ . This will be seen by the formulae given in Sec. 11. 


The relative width of the cross section 


will, how- 
max 

ever, vary considerably with @ , the value of which cannot be 

predicted with any degree of certainty. 


We have therefore chosen to use our theory to calculate 
such profiles of equilibrium as can carry the same discharge 
as in the test in an area of the same size and with the same 
Sai This is done for a number of values for i and the 


dimensions and the shape of the calculated profiles are compared 
with those of the tests. We also calculate for each value of @ 
the values of p = and thereby find k and further the value 


of the slope I, which latter quantity can be compared with the 
measured values. I is probably the least accurate of the measured 
values because uniform motion only existed over a length of about 
10 m, so that a very strict agreement between calculated and 
measured values of I can hardly be expexted. 


That value of eP is chosen which in the most essential 
features gives a good agreement in the shape of the profile, and 


4 
ta 
| 
14 
| 
| 
| 
| 
| 
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it is at the same time found that in this selection the values of 
I agree quite well. 


It also becomes possible to compare the calculated value 
of the roughness k at the bottom with the diameter of the sand 
grains, and a ratio of these quantities can be found which seems 
to agree quite well with statements from elsewhere. * 


10.3. DERIVATION OF FORMULAE SUITABLE FOR ANALYSIS 
OF MODEL TEST DATA. 


As known quantities are considered 
9,9» Tray Qand F , and we want to find 


y 
p = eo ) k, B and I, whereupon the complete cross 


section can be drawn and compared in detail with the observed 
one, and the computed value of I can be compared with the 
measured one. 


We have 
(102)  * 3 as a known quantity and start from the formulae 
in Table 8.1./15/ 
0.15324 


(91)(91c)d) = 3.1815. p whence 
max 
2. max 
6.5257 


y 
(103)]) p= = = » where 


(91e) = 1.013565 - 0.063437 log” 


From (103) p can be found; according to 8.1. (9le) Avp 
is a function of p, so we start with Avp = 1 and use successive 
approximation. p is seen to be independent of f ‘ 


When p is found, the other unknown quantities can be found 
directly: 


From Table 8.1./14/ we find: 
0.1091 
F p 


(90)(90d) 


fe A 


fp 
ymax 


*) 


For tests Nos. 1, 2 and 3 are found respectively 2.91, 1.84 
and 3.00. 


0 
| 
| A 
2.5 - 3.1815 vp 
> 
4 whence: 
t 
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(90c) A 


where 


= 1.03153 + 1.63688 cot@m - 9,02187 
+ 0.000431 cot®@ 


(90e) = 1.007083 - 0.031907 log* ( 


a is seen to vary with @, because Are is far from 
being a constant. 
For the hydraulic roughness k we have 


J 
ke mex 
P 
From Table 8.1./17/ we find 
0.0890 


(93)(93d)} B = Ang bp Ymax , where 


(93c) Abe = 1.29443 + 2.23919 cotp - 0.03069 cot*p 
+ 0.00109 cot*p 


(93e) Ay, = 1.0047812 - 0.0212087 - log*( 


From Table 8.1./18/ we find 


A t 
(94) (04a)] 1 Aj, » where 


(G4c) Aye = 1.27068 + 1.51750 cotp + 0.00519 cot*p 
- 0.00184 cot*p 


(94e) Aj, = 1.0069177 - 0.031130 log? ( 


10.4. NUMERICAL CALCULATION OF y_ .,, B, I, and k FOR 
MODEL TESTS. 


The formulae developed in Sec.10.3. are used. As units 
we use grams and centimeters. 


For all three tests is reckoned 
2 
Tmax * 9-075 kg/m* = 0.0075 gr/cm*. 
Specific weight y = 1 gr/cm®. 


max 
“VA 
fp 
: 
¥ 
| 
: 
: ; 
| 
3 
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Acceleration of gravity g = 981 cm/sec’, 


« 3.7128 cm/sec , 


y ra 6.5257 


and 10.3 (103) becomes 


1 
vp 


which formula is valid for all three tests and independent of 
the choice of 


For the calculations are used the other formulae from 
10.3 besides (103a). 


10.4.1 Test No 1. 


Measured: + = 0.0075 gr/cm%, 


Q ‘i 13.8 L/ sec = 6900 cm? /sec, 


446.2 cm* 2 
F = = 223.1 cm*, 


I = 1.35 0/00 = 0.00135, 


» = §.6 cm, 
p 100 cm . 50,0 cm. 
(102, 1) = 30.928 cm/sec. 
Ym 
575 1.4335 and from 10.4(103a) 
(103a, 1) 


where A, is taken from 8. 1(9le) 
p = 13.58 gives 
Avp = 0.96059 


p= = 13.58, final value for all 


A 85257 _ 9 7696 
vp 


a 
= = 
P = | 
| 
bs = 
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For use in 10.3(104) we find 
0.05455 0.05455 


Pp = 13.58 = 1.1533 and 
13.58 
Ay, = 1.007083 - 0.031907 - log” ( ) 
Ay, = 0.981200 ; VAtp = 0.99056, 
and 10.3(104), with F = 223.1 cm, becomes: 
223.1 
Ymax 
13.075 
max 


This expression for Foun is valid for Test No.1 and all 
values of @ . 


For use in 8.1(93) we find 
0890 _ 13, 589-9890 _ 1 2620 and 
13.5 
A, = 1.0047812 - 0.0212087 log* ( = 0.987819 
and 8.1(93) becomes 


(93, 1) B = 1.24663 Ane . (valid for all g's). 


Ymax 
For use in 8.1(94) we find A; = 0.981693. 


tT 
= 0.0075 cm 


F = 223.1 cm? 
and 8.1(94) becomes: 

fs - 0.0075 - 1.3270 - 0.981693 - 
(94,1) I = 4.37934 - 107°. 


Ymax 


For the calculation of Test No 1 for different values of p 
we now have p = 13.58 and (104,1) , (93,1) and (94,1). 

The respective values of Ay P Ang 

by 8.1(90c) , 8.1(93c) and 8.1(94c). 

The results found for m = 18°, 19° and 20° are compiled 


in Table 10.4.1/24/, where the calculated and observed values of 
» B and I can be compared, and the calculated value of k is 


given. 


and Aj ? are found 


Aip Ymax (valid for all @ ’s). 


3 
4 
4 
: 4 
3 
4 
4 
: 
a 
i 
2 
4 
4 
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Table 10.4. 1/24/ 


87 


Test No.l. p = 13.58. 
7 * 18° 19 20° 
A 2.4238] 2.3689 | 2.3179 
13.075» 
Ymax 7 = (cm) [5.39 | 5.52 | 5.64 
VAte 
| observed (cm) 5.6 
Ang 7,.9270| 7.5624 | 7,2375 
B = 1.24663: (cm) [53,31 | 52.06 | 50.90 
observed (cm) 50 
Aip 5.9458| 5.6857 | 5.4502 
-5 -5 
I = 4.37934- 10 ”- 140.5 | 137.4 |134.6 
I observed 135 
y 
k (cm) }0.397 | 0.407 [0.415 


10.4.2 Test No. 2. 


Measured: 


t = 0.0075 gr/cm? 


@ = 33.6 I/sec _ 26800 cm/sec, 


2 
a 1444.0 cm = 722.0 cm?, 


I = 0.79 o/oo = 0.00079, 


9.0 cm, 
B = 138 = 99.0 cm 
The calculations are made similarly to those applying to 
Test No.1. 
We find: 
_ 2680 
( 102, 2) Vin = 50 = 37.050 cm/sec 
1.7173 and from 10. 4(103a) 
(103a, 2) p = Fh whence 


Avp 


ll 
s). 
id 
3 
a 
ed : 
f 
Ss 
= 
§ 
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p = 34.698 
(104, 2) 22.148 tor all 9's. 


(93,2) B= 1.3709 y 


for all 


(94,2) I = 1.5244 107° 


The results are given in Table 10.4. 2/25/ 


Table 10.4. 2/25/ 
Test No, 2. p = 34.698 


Mio for all P's. 


ae 17° 18° 19° 
2.4833| 2.4238] 2.3689 
_ 22.148 “te 
(cm) | 9.14 | 9.35 
Ate 
—_— observed (cm) 9.0 
Ang 8.3283] 7.9270] 7.5654 
B+ 1.3709 (cm) [101.83] 99.30 | 96.97 
B_ observed (cm) 99 
6.2345] 5.9458] 5.6857 
-5 
[= 1,5244- Aig | 10 84.8 | 82.8 | 81.0 
1 observed 79 
y 
k = > (cm) .257 | 0.263 | 0.270 


10.4.3 Test No, 3. 


Measured: + = 0.0075 gr/cm’, 


q = 87.0 I/sec 43500 cm3/sec, 


2 
F = 1237.5 cm?, 


I = 0.65 o/oo = 0.00065, 
= 10.5 cm, 


Ymax 


| 
| | 
| 
i 
: 
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B- = 140 cm. 


The calculations are made similarly to those applying to 
Test No, 1. 


We find: 
43500 
(102, 3) Vm = {337.5 = 35.152 cm/sec 
375 1.6293 and from 10.4(103a) 
24.180 whence 
(103a,3) p= 3557 
A 
vp 
p = 26.10 
29.517 
(104, 3) for all 
f 
? 
(93, 3) B = 1.3332 - Ane for all@ ’s 
-§ 
(94, 3) I = 0.85874 . 10°. Aig’ for all P ’s. 
The results are given in Table 10.4.3/26/ 
Table 10.4, 3/26/ 
Test No.3. p = 26,10 
14° 15° 16° 
Ar 2.6969; 2.6189) 2.5481 
29.517¥ 
a (cm) 10.94 | 11.27 | 11.58 
A 
observed (cm) 10.5 
Ane 9.8520] 9.2804 | 8.7764 
B = 1.3332 - Ane (cm) 143.76) 139.45 | 135.04 
B_ observed (cm) 140° 
Aip 7.3310) 6.9199 | 6.5571 
-5 -5 
= 0,85874- 10° . Aip 68.9 66.9 65.2 
observed 10 65° 
y 
k = (cm) |0.420 | 0.431 | 0.444 


| 
4 
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10.5. ADJUSTMENT OF THEORETICAL EQUILIBRIUM PROFILES 
TO THOSE OF MODEL TESTS. 


From Tables 10.4.1/24/, 10.4.2/25/ and 10.4.3/26/ we 
already get a rough idea of the value of @ that will give the 
best adjustment of the theoretical equilibrium profiles to the 
observed ones. 


To compare the profiles more closely the theoretical 
profiles are drawn according to the co-ordinates given in 
Table 8. 2/21/. 


The observed profiles, the co-ordinates of which 
are given in Table 10.1/23/, are then reduced to the same scale 
as the theoretical ones, and superimposed correctly upon the 
theoretical profiles, i.e. with coinciding axes and water surfaces. 


This has been done for three different values of Pp for 
each tested profile. 


That value of @ is chosen which in general gives the 
best conformity between calculated and observed profile for 
each test. 


The following values of @ were chosen: 


Test 1: = 19° 
- 83: 
@ = 14°. 


The results are shown in Figs. 8, 9 and 10. 


It is seen from the Tables /24/, /25/ and /26/ that the 
agreement of I corresponding to this choice is quite good 
considering the uncertainty of these observations, as mentioned 
in SEC.10. 2. 


As seen in the figures the agreement is considerably 
better for Bests 1 and 2 than for Test 3. 


If we study the data for Test 3 we find, however, that its 
mean velocity Val is smaller than vi for Test 2 although or 


is ~ 17% greater in the former case. This is in disagreement 
with the well-known fact that for equal tractive force a greater 
mean velocity is needed to move the sand grains when the depth 
is greater. It therefore seems as if some uncertainty must exist 
in some of the observed data. 


It does appear from all three tests that the angle of 
internal friction @ , which must be introduced in the theory put 
forward to get a tolerably good agreement, is considerably 
smaller than would generally be expected. 


If we remember, however, the fact that the sand grains 
in a state of incipient motion hardly exert any normal pressure 


ee 

| 
3 
ex 
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on the underlying ones and that the angle @ is known to de- 
crease with decreasing density of the sand, it is hardly sur- 
prising that @ is found to be small. 


In Figs. 8, 9 and 10 the slope of the banks in the model 
tests is seen everywhere to be greater than for the theoretical 
profiles, which fact might seem to contradict the reasoning above. 


It does seem possible, however, that the capillary tension 
in the sand above the water surface may be able to give the 
banks a certain cohesion causing them to be locally steeper than 
in the case of a true cohesionless material. 


A similar effect in a full-size channel or natural water- 
course would of course be relatively insignificant. 


We, therefore, provisionally conclude that the angle of 
internal friction @ for a sand material like that used in the 
model tests should be estimated at 18 to 20°, but it is clear that 
a further study of this question is highly desirable (see Sec.14 ). 


11. FORMULAE FOR THE DIMENSIONS OF 
EQUILIBRIUM PROFILES FOR GIVEN 


In the formula 8.1(89) for Q’ in Table 8 1/13/ we insert 


A 3.1815 - A 
8. 1(89c) ae 

and 
8. 1(89d) Bap 0.26234 = constant 

and get 

2. 26234 
Aap 

where 

8.1(89e) AQ, = 1.020619 - 0.095322 log* 


and where Ate is found from 8.1(90c), Table 8.1/14/. 
Solving with respect to p we get 


p. 1,009114- [1 - 0.041283 log? ( 
Q 


0. 44202 


~ 
2.5 — ce. 
where 7 


8. 1(90c) Atg = 1.03153 + 1.63688 m - 0,02187 m? + 0.000437 m® 


| 
+ 
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and 
m= cot@. 


In the first approximation the factor next to p on the left 


side of (106), viz: 1.009114-[1 - 0.041283 log” ( rrP-aem)] is 


neglected, whereby a nearly correct value of p is found; p is then 
found correctly by successive approximation. 


For 20 £ p < 600 the said factor will hardly deviate more 
than 1.25% from unity. 


y 
When p = yo is found we get directly 


(107) k p. 


For the half width of the water surface, B, we get from 
Table 8.1/17/, formula 8.1(93) with 8. 1(93d) 


Bip = 0.0890 = constant 


Bek: 0890 


Ab, 
where 
8.1(93c) Ang = 1, 29443 + 2. 23919 m - 0.03069 m*+ 0.00109 
and 
m = cot 
further 


8.1(93e) = 1.0047812 - 0. 0212087 ( - 


As p is known from (106) B can be directly found from 
(108). 


_ From Table 8.1/14/, formula 8.1(90) we get with 


Br = 0.1091 = constant for the area of the half cross section 


2.1091 
where 


8.1(90e) Ag, = 1.007083 [i - 0,031682 ( )}. 


p is known from (106) and F can be directly found from 
(109). 


We want, however, to find a more explicit form of the 
formula for F and therefore insert the expression (106) for p 
in (109). 

We find : 


: 
4 
(109) 
‘ 
= 
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=, 


2.1091 
9114+ [1 - 0.041283 ( 


1.019223-[1 - 0.087070 log* ( 


By using the formulae 6.4(77 through 80a) we further find: 
1.019223 - 0.087070 log* ( 


1.007083 fi - 0.031682 
1.012051°[1 - 0,055388 log” 


and (109) takes the form : 
13547 


0.067735 
Ar 
max 
2. 
= 
2.97877: [I - 0.055388 log? ( 
where 


8.1(90c) Ate 1.03153 + 1.63688 m - 0.02187 m* + 0.000473 m® 


nu? 


and 
m = cot@. 
0.067735 
In (109a) the factor Ate takes the following values 
for the four values of P : 
0.067735 
At@ 
15° | 1.1393 
20° | 1.1206 
25° | 1.1063 
30° | 1.0947, 
and it is possible to put with good approximation 
0.067735 
Ay - 1.0560: [1 + 0.02178 cot@], 


where upon (109a) takes the following final form : ohne 
0, 
0.13547, 
F -0:354509 k [i + 0.02178 cotg] 


1 - 0.055388 ( ) na 
From (109b) it is clearly seen that F is practically 


independent of p = ae for actual channels, where p will 


0.932265 
4 
| 
4 
38 
a 
i 
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probabl y be of the order of 100 and that it varies only slowly 
with @. 


For practical calculations it should consequently be 
advantageous to judge roughly the value of p from (106) and then 
use (109b) to find F. 


From (106) it is seen that p varies approximately pro- 


portionally to "0 88404 ang thereby it is seen from (107) that 
is approximately proportional to 11506 whereas (108) and 


(106) show that B is approximately proportional to 


k k 93728 
7,0. 00404,1.0890 96272 
) k 
The dimensions proper of the equilibrium profile 
consequently are only slightly influenced by an uncertainty of k. 


11.1. CHOICE OF THE ANGLE OF INTERNAL FRICTION @ AND 
THE EQUIVALENT SAND ROUGHNESS k. 


11,1.1.9. 


As to the choice of @ we refer to the remarks made in 
Sec.10.5 . 

The value of ® must probably be lower than those values 
used in calculations of soil mechanics, and the results of the Vienna 
tests make it probable -that @ = 18° to 20° will be appropriate. 
Further studies of this question do, however, seem necessary, and 
we refer to the proposals made in the conclusion, Sec. 14. 


hk. 
As to the value of k, the equivalent sand roughness, which 


must be introduced in the formulae above, several view points can 
be advanced. 


coefficient M in the formula v_ = R2/3. 


If an empirical knowledge of M for channels of a similar 
character is available, it is easily proved by comparing the Manning 
formula with the results of Nikuradse for the resistance number f, 


that 
25. 84 


/sec) 


or M 6 
(110) k = ( 35-37 ) (k in meters). 


11.1,.2.1. Determination of k means of the 

(metric units). 


ind 


1D 


nna 
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If Kutter’ s n is known we have 
(111) M =i. 


The above-mentioned results of Nikuradse concerning the 
resistance number f agree perfectly with the results obtained by 
pure calculation on the basis of the logarithmic velocity distribut- 
ion, as it is done in this paper. 


From the logarithmic velocity distribution in a definite 
normal to the bottom of total length Zo it is easily found that the 


velocities YD and Yq at points p-z, and q: Zo from the bottom will 
have the mean value Ven equal to the mean velocity in the normal 
concerned, provided that 


(112) p-qs + = 0.1352 (e basis of natural logarithms). 
e 


v.79 
(113) vn 


The hydraulic roughness of the bottom in the neighborhood 
of the foot of the normal is found to be 


29.7 - Zp 
(114) as ( p>q). 
Vv. 
p “q 


The values p = 0.90 and q = 0.15 e.g. satisfy (112), and 
if the velocities are measured at distances 0.90- Zo and 0.15 Zo 
from the bottom we get from (113) and (114) 

26.7 - 20 


(114a) k = for p = 0.90and q = 0.15. 


6 


By taking normals distributed over the whole width of the 
channel we can get an idea of the variation of k transversally to 
the axis of the channel and obtain a mean value of k. 


If empirical values of k inthis way are known from 
channels of varying character, a sufficiently good estimate of k 
for planing purposes can probably be made. 


ii,1,2,2,_k found from observation of velocity distribution _ 
in_a_normal to_the bottom. 

| 
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an 
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Provided that the channel considered is shaped as an equi- 
librium profile a relation can be shown to exist between 
Ymax vm 
p= and as well as between p and , the ratio of 
max 
mean to maximum velocities for the total cross section. 


If the cross section is measured in detail 9€ - 


can easily be found, and the corresponding value of p = a 
can be read in a diagram; we then find 
Ymax 
Pp 
A possible knowledge of 


Vmax 


can similarly be used for 
finding k. 
This method is presented in detail in Sec. 12,(see Fig.11). 


11,1,2.4. Assuming the same value of k as emRcahly 


found for natural watercourses, viz. = » where 


This method will be developed in detail in the following 
Sec.11.2, where formulae for p, \ Pe B and F are presented 


that do not directly contain the roughness k. 


11,2 DETERMINATION OF DIMENSIONS FOR GIVEN VALUES OF 
o te 9 AND @ AND FOR k ASSUMED TO CORRES- 


POND TO ROUGHNESS OF NATURAL WATERCOURSES, 


As it will perhaps sometimes be doubtful which value must 
be assigned to k in the formulae of Sec.11, we think it useful to 
develop formulae not containing any value of k, which can be 
chosen more or less arbitrarily. 


From a study of many hundreds of observations of 
velocities in natural watercourses together with their hydraulic 
radii R and their slopes, we have found statistically that k must 
be independent of the slope I and is a function of R alone. 


11.1,2,3,.__k found by means of degree of fullness, 
Vv 

m 

or 

max 
= 0.425 meters or 1.395 feet. 


We have found 


t= 
where j 

e* 0.425 meters or 1.395 feet. 


For small values of R, (R<1.0 m), e.g. for mountain 
torrents, the dispersion of the observations in relation to the law 
(115) is considerable, Sut for larger rivers in plain country the 
agreement between the observations and (115) i8 very good. 


It must naturally be remembered that k itself varies during 
floods and that the same value of Q is known to correspond to 
different depths by rising and falling waters. No great exactitude 
can consequently be expected. 


When we set up a velocity formula based on (115) and on 
the ordinary assumption that shearing stresses are uniformly 
distributed along the perimeter we find 


R 
(116) 8,83+[0,750 + log 


This result is in good accordance with some of the better 
formulae for natural watercourses, e.g. those of Hermanek. 


This question is further treated in Sec. 12. 2. 


In the following we therefore use the assumption (115) and 
(115a). 


(115) is written as follows: 
1 56 


7, \ 2.785 


From Table 8.1/16/ we find 
0.0210 


max 


_ ¥max 
by multiplication of 8.1(92) by p = «CWE get: 


(117) App: Py where 
8.1(92c) Are = 0.732 and 
8.1(92e) AL, = 1.000 
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From (115b) and (117) we get 


c 
k 0.358974 —(0.3665125 
A,») P 


y 
By division of p = ae by (118) we get 


Ymax -0.358974 _0.6334875 

Using (118) and (119) k and Ymax @Pe now expressed by p. 
In Sec. 11 we found : 


11.(106) 1.009114 - [1 - 0.041283 log” ( 


1 Q 0.44202 
0.88404 
max 
2.5 -\/——— - 3.1815 : A 
9 


(118) is raised to the power 0.88404 and multiplied by 
11,(106) and, when solving with respect to p, results in: 


§53887 


479316 


0.717948 


max 2 
2.5 = 3.1815 Ae | 


Introducing the values 8.1(92c) and 8.1(92e) 
Ale" 0.732 and Arp = 1,000 we get: 


p [1 - 0.06107 log” | 


] 0. 653887 
= 0.39987 


Q 
2.5 Aral 


In (120a) p is expressed by known quantities. The paren- 
thesis on the right side is dimensionless. Note that the variation 


of Ate with ed is considerable. 
For rough calculations we take for instance 


Q: (Ag: A, 
| | 
| 


6 
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and estimate p = 110; i.e. log* - 0. 


We then get: 


0.6539 
Pp = 0.1333. (P = 20°) 


max 2 
2.5 9 Cy 
valid for every consistent system of units. 
In the Metric System we have 


T3078 0.3527 
0.425 °° 
0.6539 
p = 0.408 (@ = 20°) 
2. 5\|_max 
3 
In the FPS-System we have 
1 1 
0.6539 
p = 0,0863-| (p = 20°) 
2. 5: max 
3 
From (118) we get 
-0.358974 —-0.3665125 


The equation (120) for p is raised to the power -0.3665125 
and multiplied by (118a), whereby we get 


‘ h. fi - 0.041283 - log* ( 


k Ago) | Q 239658 


542189 


(121) 


p. 
| 
| 


100 ; SEC. 11.2. 


In formula 11.(109a) we found ; 


yO: 13547 0, 932265 


Ay 


? 2.5. 


2.97877 [1 - 0.055388 


When (121) above is raised to the power 0.13547 and 
multiplied by 11.(109a) and observing that 
2 - 0.13547 = 1.86453 = 2- 0,932265 we get 


11.(109a) F = 


2.899799 
| Q 


A 


73450 


(122) 


2.97877 | — 0.055388 log "GA 55 | 


We put: 
0.073450 


and find by formulae 6.4(77 through 80a) 


1 - 0, 055388 log’ = 1 - 0.052356 log* (rrf ) 
1 - 0,00303224 


which result is inserted in the denominator of (122). 
We further have 
0.732 , A. 1.000 
0.071940 


= 0.97780 
and find the new form of (122) 
A 0.100201. Q 0.899799 
4 
2 


c,? 2.97877 - 0.97780-[1 - 0.052356 log’ 


4 

i 
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0.100201 


The variation of the factor Are with@ is only 


slight, and we get with good accuracy : 
(123) 190208 0956 (cotp)”: 
With (123) we finally get from (122a) 


0 899799 
0. . 2. 
2.5V . c,? 


- 


- 0.052356 - log* ( ) 


It is clearly seen that the variation of F with @ is only 
slight. 


For rougher calculations we put 


0766 ~ 
0.0766 - 


p= 20° , (cotp 1.0812 and 


p ~ 112° t.e. log* ) = 0 


and get 


(p= 20°) 


0.9 


(Pp = 20°) 


This expression can be used for every consistent system 
of units. 


In the Metric System 
= 0.425m = 0.8427 4) 
0.9 


2. —— 


In the FPS-System 


* 1.395 feet , = 1,069 (feet? ? 


F 0,435. (Pp = 20°) 


« 0.9 
<= 0.407: 
- c,* 
k 
~ 0.2 
FF 0.407 - 
ij max 
5 2.5 
max 
| 
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From (119) we get with 


Are = 0.732 Arp = 1,000 
y 
mex 1.11852 633488 
k 


When p is found by one of the formulae (120a)b)c)d) 


is directly found by (124) . 
From Sec.1l. we have : 


11.(108) B= Aye ik 
We find by division by (118) 


B. 0.72249 


k 


The variation of B with @ is considerable so that the 
correct values of Ang should be taken from Table 8.1./17/, 


formula 8.1(93c). 


The variation of A, will probably only be small 
(cf. formula 8.1(93c)). P 


By the successive use of the formulae 


y 
(120a)b)c)d) for p= and 


bp ° 


(122a)b)c)d)e) for F and 
(124) for and 
(125) for B 


all the principal dimensions of the equilibrium profile can be found 


and the complete profile can be drawn using the tables of the 
co-ordinates ) and § , Table 8. 2/21/. 


It must not be forgotten that Q, F and B refer to the 
half total cross section. 


| 
{ 
| 
é 
— 
. 
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11.3. ILLUSTRATIVE EXAMPLE OF CALCULATION OF DIMENSIONS 
The Metric System is used, 


Given : 


Q= m®/sec = 400 m*/sec for half cross section, 


Tmax 0-15 kg/m*, 
1026 kg/m’, 

Pp = 20° 

and 

k ~ as for natural watercourses. 


For use in the formulae we calculate : 


t 
max mae = 0.001435 (m/sec)’, 


0.03785 m/sec, 


. 0.4225"°9 = 3980 - 0 4606 = 1833 


and further have 


11, 2(115a) c, = 0.425 m, 
11, 2(122d) F = 0.343 - 1833 = 628 m? 
Total area: ZF = 1256 me 
0.6539 
107-8156 9 49959-6539 _ 419. 9 5700 
2.5 = 234.5 


y 
11.2(120c) p = = 0.408 - 234.5 = 95.7 


p?- 6335 95 70-6335 _ 17 95 


11.2(124) 0.425 1.1185 - 17.95 = 8.53 


| 
| 
t 
max 
2.5° 
0.9 3.6 
max 
5) 
} 
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= 0.0891 m, 


P 


7225 0.7225 


95.7 = 27.00. 
From Table 8.1/17/ we get 


Ang 7.2375, and 8.1(93e) A, = 1.004781 - 0.02121 log* 


106.46 
1.1125 ; log1.1125 = 0.0462 ; log® 1.1125 = 0.002135. 


Ay, 1.004781 - 0.000045 = 1.004736 


11, 2(125) B = 0.425 1.1185 7.2375 27.0+° 1.0047 = 93.2 m 
Width of water surface of total section : 


12, FORMULAE FOR THE MEAN VELOCITY v_ IN 
EQUILIBRIUM PROFILES. 


12.1 Sans AS A FUNCTION OF HYDRAULIC RADIUS R. FIXED 


VALUE OF k. 
From Table 8.1/16/ we have 
8.1(92)c)d) R’ = = 0.732 210 


8. 1(95)c)d) 


or 
(127) 


From Table 8.1/15/ we have 


Vv 
8. 1(91b)c)d) v’ = m = 9; 1815 15324 A 


where g 


8.1(91e) = 1.0135645 - 0. 0634373 ( ) - 


R 0.732 pi: 210 

1 

or 

R 210 

: From Table 8.1/19/ we have 

max 

1-160: 
max 
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0.01 0.15324 
Vin 2:5 3.1815 + 1,160 + p “Wert p 
or 


0.16324 
(128) vi. = 2.5 + 3.1815 1.160 p RI 
For the calculation of Avy we find by means of (126) 


1 
log * 10K | 


1 


from which 


and 


0.0634373 - log ( = 0.043329 ( » 
which inserted in 8.1(9le) gives 


5) 
35. 


: R 
(129) 1.013565 - 0.043329 log? ( 
We further have from (126): 
| 0.16324 
| (130) 0.1694 
0.732 k 0.7320: 14008 
Inserting (130) in (128) we find 
m 07390: 134909 k 
where Av is taken from (129) and 
0.7329: 134909 = 9 
We finally get 
(13a) v 
0.134909 
= 9.6228 - [1.013565 - 0.043329 log 


| or 


© 
f 
3 
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k 


(13la)b) is valid for every consistent system of units. 
In the Metric System we have 


Ve = = 3.1321 , 


which gives 

30.14 0.6349 R | 

(131c) R I *11.0136 - 0.0433 - 
m k . [ 


Herein the last parenthesis (i.e. A, ) is practically 
always near unity. P 


12.2 an AS A FUNCTION OF HYDRAULIC RADIUS R, ASSUMING 
k AS FOR NATURAL W 


From 11, 2(115)a), valid for natural watercourses, we have 
= 
where 
11,2(115a) c, = 0.425 meters or 1.395 feet. 
This is inserted in 12.1(13la) and we get 
R 0. 21046 
(132) 9.6228 - Avp 


where 


1.56 
(132a) = 1.013565 - 0.043329 


We have 
1.56 


which gives 
R 
(132b), AY, = 1.013565 - 0.105445 log*® 


(131b 
1 

108) = log} ( = 1.56 log ( 

lo | 
: | 
| 
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Inserting (132b) in (132) we finally get 
R 0, 21046 


(183) v,, * 9.6228(8.) Vert 105445 lof 


or 

"ie 
9.6228 Vg 0.71046 [ R 
|1,013565 - 0, 105445 log* 


which formulae are valid for an arbitrary consistent system of units. 
In the Metric System we have 
Ve = = 3,1321 m!/2/ge¢ 
0.425 m ; 31.7- 13.48, 


0.4259: 21046 0.8350 


and (133a) becomes 
(133b) 36.0953 - 71046 91357-0. 105445 


This formula has much resemblance to the empirical 
formulae of J. Hermanek. 


Hermanek gives for R = mean depth of section and 
1.5m<R<6.0 m the following formula 
0.75 1/2 
(134) 34- R - I 


The last factor in (133b) will be near unity for such 
values of R. 


12.3 Vn AS A FUNCTION OF y FIXED VALUE OF k., 


If we want to express the mean velocity by Yuna instead 
of by R we insert R from 8.1(92) in 12.1(128) and get 


0.5 _0.17374 
v. = 3.1815 - 1.160 0.732 p y 


(135) v_ = 7.89379 - -I° A 


where 


y 
(135a) AY, = 1.013565 - 0.063437 - log* 


or 


‘ 
| 
| 
| 
. 
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7.89379 Ve 0.67374 1/2 
vm *[ | Avy 


k 


(135)a)b) are valid for an arbitrary consistent system of 
units. 


In the Metric System we have 


Va 3.1321 m!/?/sec 
and 


24.7241 0.67374 1/2 
Avp will generally be very near unity and the similarity 
between (135c) and the Manning formula 


becomes striking; it must, however, be remembered that (135c) 
contains » ae and not R. 


13. THE RATIO OF MEAN VELOCITY Ym TO MAXIMUM 


VELOCITY ax: 


m max 
13.1 AS A FUNCTION OF p = 


max 
From Table 8.1/15/ we have 
Vv 
8. 1(91b)c)d) 3.1618 . po 
max 
max 

where 

8. 1(91e) A 


vp = 1.013565 - 0.063437 log* ) 


For the maximum velocity Vsnums which occurs at the surface 
where the depth is Fh... have directly 


v y 

(136) mex 3.392 +n - 3, 392 [1 + 0.67883 log p| 
25 max 
We VR 8.1(91b) by (136) and get 


0.15324 
Vm 0.950665 p - 0.062588 log* 
Vv 7 I + 0.67883 log p 


} 
3 
i 
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The ratio of mean to maximum velocities is 


max 

consequently a function of p alone and quite independent of P . 
The denominator in (137) 

v 


1 + 0.67883 log p = 


2.5 max 


can be developed in a form similar to that of the numerator in 
(137). 

The calculation is made in the same way as that used 
for the functions A. in Sec. 6.3, 


We find as a first approximation 
(138a) 1 + 0.67883 log p = 1.31155 p’' 17595 + 9.925 % 

for 24.25 p < 500. 

The logarithmic deviations ep from the simple power 
formula are then treated as indicated in Sec. 6.3 and we find 
(138b) 1 + 0.67883 log p 
= 1.323712 19939. 9, 041934 log? |? 0.104 % 

for ~ 18<p<~ 950. 

Inserting (138b) in (137) we get 


0.02789 
Vv 


max (1 - 0.041934 log* 


The parentheses are divided by means of the formulae 
6.4(77 through 80a) and we finally get the following form of (137a) 


Vv 


max 


13.2 DEGREE OF FULLNES # OF THE CROSS SECTION OF 
EQUILIBRIUM PROFILES, 


We define 3€ as follows 


F 

F Vmax 


max 


p 
| 
| 
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From Tables 8,1/14/ and 8,1/17/ we get 
70.0890 
(93b)d) B’ = p 
and by division 
A A 
F’ f 0.0201 f 
B’ bp 
Are 


For the ratio 7—— we get from Tables 8,1/14/17/ the 
following values oP 


A |0.739031 |0,.742356 |0.745156 | 0.747759 
be 


The variation with p is so small that we take as an 
average A 


(139a) mS = 0.7424 = constant , nearly corresponding to 
Pp = 20°. 
For Ay, and Ay, we get from Tables 8.1/14/17/ 


(90e) Ap, = 1.007083-[1 - 0. 0316826 log* ( rag 
(93e) Ay, = 1.004781-[1 - 0.021078 log* | 


The ratio of the two parentheses is found by means of the 


formulae 6.4(77 through 80a). We finally get 


(139b) = 1.002290: - 0.0105748 log | 


and by inserting (139a) and (139b) (139) takes the faqllowing form 


2 = 0.7441 - [1 - 0.010575 log” 


| | 
: | 
| 
. 
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Vv 
13,3 RELATION BETWEEN AND 2€. 
max 
We found 
max 


13.2 (139¢) = 0.7441 0, 010875 10¢" 
Vv 


These equations represent 


= and @€ as functions of 
max 


Y max 


p = ~~ — ; the relationship between and @@ is given by 


max 
Fig. 11 with the parameter p on the curve and the figure can there- 
fore also be used for the determination of p and thereby of k when 
one of the other two variables is known, 


In order to eliminate p between 13.1(137b) and 13, 2(139c) 
by raising the latter to power 1.38756 and dividing we find 


Vv 
) 1 - 0.020654 log* ( any ) 


Ox = 1.082324 
(140) 1 - 0.010575 log* ( ) 


In (140) the fraction on the right side is treated according 
to the formulae of Sec.6.4 and we get 


Vv 

(140a) = 1,082352-[1 - 0.005981 log* ( )| 
To eliminate p in (140a) we use 


13.2(139¢) = 0.7441 - - 0.010575 log” )] - 


Taking logarithms in 13.2(139c) and noting that the 
parenthesis does not deviate much from unity we find 


(141) flog rrP yas H 0. 22851 log ( 


49.751 log (1.22252). 
We substitute 


(142) € = 1,22250€ - 1 


and noting that € will always be small compared to unity, (141) 
becomes 


fice rrr - 0. 22851 low - 21.606 -E, 
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Solving this quadratic equation we find 
log ( ) 2.1861 [1 - 19. 4706 | or 


(141a) log ( = 2.1096 [1 - 1.0372 - 19, 4796 |. 


(14la) is inserted in (140a) and we finally find 


v 
1.02285 0. 5904E+0, 05845\i-19 4706 | 


max 


where € is taken from (142). 
(140b) gives directly the relation between @€ and — 


We find for instance ¥ Vmax 
2e m 
“max. 
0.77 0, 7478 
0.80 0.7929 
0. 84 0. 8479 


It must be remembered that (140b) is exclusively valid for 
equilibrium profiles. For other channels, e.g. with fixed bottoms, 


max vm 
no definite relation between @€ and p = “Ee | COR exist and 


max 
can vary within considerably wider limits. 


For an equilibrium profile it will consequently be possible 
to determine the mean velocity by measuring the maximum 
velocity and sounding the complete cross section of the channel, 
i.e. determining #. 


This should be an advantage in tidal channels, where the 
discharge cannot be considered steady for a time sufficient to 
carry out velocity measurements at a sufficient number of points 
of the cross section. 


2 depends not alone on the area of the section and its 
maximum depth which will probably both be well defined, but also 
on the width of the water surface. 


This width will e.g. for natural streams often be influenced 
by vegetation on the banks and thus be narrower than for the 
equilibrium profile, even if the cross section as a whole is near 
the equilibrium profile. The area can then be considered correct 
and the width of that equilibrium profile which as an average 
gives the best adjustment to the profile in question can be used by 
the determination of 9€. This is what has been done in Sec.10 by 
the study of the Vienna model tests. 


of 

| 
| 
| 
| 
| 
| 
| 
| 
| 
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SEC. 14. CONCLUSIONS AND SUGGESTIONS 
FOR FURTHER STUDIES. 


The topics studied in this paper would be of interest for 
the problem of stable channels in general, but they do not cover 
the whole problem. 


The shape of the profile which is stable against erosion 
is found, but only provided that no transport of bed load occurs; 
neither has the question of silting up of the channel been taken 
into account, 


It is hoped however that the present study will be useful 
if the problem in all its aspects is to be solved. 


The whole study is carried through theoretically on certain 
assumptions which are clearly stated, and further work in this 
field should try better to confirm or to correct these assumptions 


by direct control or indirectly by comparing the results with 
actual observations. 


The most important assumptions made are: 
1) Distribution of t according to formula 1.(1) or 4. 2(1a). 


2) Hydrodynamic lift L proportional to t, 2.1(3) and 2.2(9) based 
on the work of Einstein [1]. 


3) Logarithmic velocity distribution as indicated in Sec.1; mean 
velocity in an infinitesimal wedge-shaped element of area 
assumed equal to mean velocity in normal to bottom at the base 
of elements, Sec.4,2. 


4) Uniform distribution of roughness k, Sec. 4.2. 
5) Constant value of angle of internal friction@ along perimeter. 
6) Principle of minimum area of cross section, Sec. 4.4. 


7) For some of the formulae for the dimensions of the channel 
and for the mean velocity an empirical relation between 
hydraulic roughness k and nydFaallé radius R, as found for 
natural watercourses, has been introduced. Thus formulae not 
containing k have been established. 


The igovels studied in Sec. 9 and the example presented 
in Fig.7 seem to some degree to confirm assumption 1) concerning 
the distribution of +. 


That this choice is better than the usual assumption 
taking * proportional to the depth, seems evident, but the question 
could be studied further by other methods, and the result possibly 
be improved; it seems probable however that such further studies 
will be considerably facilitated when an approximately correct 
solution is at hand. 
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Assumption 2) concerning hydrodynamic lift should also 
receive closer study and the size of the proportionality factor 
c «4.9 might probably be improved. It seems clear however that 
this assumption is better than the complete disregard of the hydro- 
dynamic lift, and the large size of c may perhaps explain why the 
angle@ is found to be so small as is the case (see below), 


The logarithmic velocity distribution accepted is probably 
so well established that it does not need further comments, but the 
hydraulic roughness in real profiles must obviously be expected to 
vary along the bottom; the k introduced must therefore be taken as 
an average value. 


As to the principle of minimum of area it will hardly be 
possible to give any strict proof of this assumption, but it almost 
seems selfevident and corresponds completely to many similar 
cases in hydraulics. A comparison of results with observations 
from nature seems desirable. 


The values of @ have been determined for three model 
tests carried out in Vienna,1916, giving = 19°, 17° and 14° as 
the most probable values for these tests. The values are even 
lower than the expected value of about 20°, and it seems highly 
desirable that further observations should be made. Time has not 
allowed independent model tests to be carried out, but this ought 
to be done. Further it seems probable that more insight into the 
size of @ might be obtained by three-axial tests with sand where 
the hydrodynamic lift is replaced by an upward stream of water in 
the cylinder. By such tests the effective grain to grain normal 
stress could be reduced even to zero. This can possibly also be 
effected by other means, e.g. using a liquid with a specific weight 
greater than that of water. 


A correct value of P can thus be found for effective 
stresses smaller than those corresponding to the weight of the 
grains proper. 


The fundamental results of the study are given in the 
tables in Sec.8.1 All quantities here are dimensionless and we 
have succeeded in expressing Q’, F’, v,.’, R’, B’, I’ and v,’ 


by the product of two functions, one depending on @ alone and 


y 
another depending on p = —— 


alone. 

The limits of validity and the corresponding maximum 
deviation of these formulae are given. For many practical purposes 
the simple power formulae without the correcting factors 


should probably suffice; but for the application of the formulae for 
model tests, e.g. where p = 7max/k is small, the complete 
formulae including Ay should be used, 


: 
i 
= 
ave 
| 
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For all derivations in this paper the number of digits in 
the constants of the formulae might seem exaggerated for practical 
purposes; wane however the great number of numerical 
operations that have been needed in deriving the formulae, it has 
been judged appropriate to retain a surplus number of digits in all 
derivations, and the values of the constants can then finally be 
rounded off, 


The form of the formulae is such that all practical 
calculations can easily be carried out by means of a log log slide 
rule, 


The result found in Sec.11, viz. that the area F of the 
half cross section of the equilibrium profile varies only slightly 
with the angle of friction @ and with the relative roughness, but 
mainly depends solely on the discharge Q, the limiting tractive 
force +t and the hydraulic roughness k, might facilitate the 
comparfBor of results with observations in nature. 


Where the flow inthe channel is caused by tides, the use 
of the results in Sec.13 might facilitate the observation of Q. 


Until more definite knowledge about the roughness k is 
acquired, the results in Sec.11.2, where k is taken to be equal 
to the empirically determined roughness of natural watercourses, 
might be of use. The results thus obtained should, as far as 
possible, be compared with observations in nature. 


The formulae for the mean velocity found in Sec.13 also 
make it possible to find k from available observations. 
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